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ABSTRACT 

This  study  uses  a  numerical  model  to  investigate  the  effect 
of  the  air-sea  fluxes  of  sensible  heat  and  moisture  in  extra- 
tropical  cyclogenesis  over  the  open  ocean.   The  model  is  a 
sectorized  version  of  the  UCLA  general  circulation  model  in- 
cluding the  Aratewa-Schubert  cumulus  parameterization  scheme 
and  the  Randall  version  of  the  Deardorff  planetary  boundary 
layer  parameterization.   Idealized  initial  conditions  are  speci- 
fied in  the  atmosphere  and  ocean  that  are  typical  of  open-ocean 
fall  and  spring  conditions. 

Adiabatic  and  diabatic  model  results  are  compared  over  a 
15-day  integration  period.   Diabatic  processes,  including  the 
surface  fluxes,  cause  a  large  reduction  in  low-level  static 
stability  during  the  initial  cyclone  growth  period.   This  re- 
duction in  static  stability,  as  well  as  the  latent  heat  release, 
leads  to  the  rapid  growth  of  wave  numbers  12  and  18  in  the  dia- 
batic model  experiment,  while  only  wave  number  6  is  present  in 
the  adiabatic  experiment.   The  growth  of  the  cyclones  is  much 
more  rapid  in  the  diabatic  experiment.   However,  cyclones  in 
the  adiabatic  experiment  attain  similar  maximum  intensities  as 
in  the  diabatic  experiment,  and  undergo  an  analogous  decay  period 
After  development  of  the  initial  cyclone,  a  variety  of  secondary 
low  developments  occur  in  the  diabatic  experiments  that  are 
similar  to  those  observed  over  the  open  ocean. 


The  role  of  the  air-sea  fluxes  is  studied  by  selectively 
removing  the  surface  fluxes  at  various  stages  of  the  cyclone 
development.   Removal  of  the  air-sea  fluxes  produces  signifi- 
cant changes  in  the  evolution  of  the  diabatic  model  cyclones. 
Removal  of  the  surface  sensible  heat  flux  results  in  an  in- 
tensified low-level  temperature  gradient  in  the  developing 
cyclone.   The  enhanced  low-level  temperature  advection  in- 
creases large-scale  lifting  and  therefore  latent  heat  release, 
resulting  in  an  intensification  of  the  cyclone.   Removal  of 
the  surface  moisture  flux  prevents  the  development  of  second- 
ary lows  which  might  otherwise  eventually  replace  the  primary 
cyclone,  as  occurs  in  the  complete  diabatic  model.   These 
model  results  are  not  verifiable  with  observations,  however, 
they  show  a  mechanism  by  which  some  disturbances  over  the 
ocean  could  develop  while  other  seemingly  similar  disturbances 
would  not. 
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I.   INTRODUCTION 

Extratropical  cyclones  are  responsible  for  the  meridional 
transport  of  heat  and  momentum  as  well  as  for  most  signifi- 
cant weather  in  the  mid-latitudes.   Simpson  (1969)  estimates 
that  7  5%  to  85%  of  the  mid-latitude  transfer  of  heat,  momentum 
and  moisture  through  the  ocean  surface  occurs  during  the  pas- 
sage of  extratropical  cyclones.   Understanding  the  general 
circulation  of  the  atmosphere  and  oceans,  therefore,  requires 
an  accurate  knowledge  of  these  large-scale  eddies  and  how 
they  utilize  the  large  amounts  of  energy  transferred  to  them 
from  the  world  oceans.   These  storms  also  cause  considerable 
destruction  (Gyakum,  1980;  Rasmussen,  1979)  and  disruption  of 
commerce.   Current  numerical  models  frequently  fail  to  pre- 
dict these  systems  accurately  beyond  one  to  two  days,  possibly 
due  to  the  inadequate  parameterization  of  the  diabatic  proc- 
esses (Sanders  and  Gyakum,  1980)  .   The  success  of  medium  and 
extended  range  forecasting  may  depend  on  the  proper  under- 
standing and  modelling  of  the  role  of  the  air-sea  fluxes  in 
the  development  of  the  extratropical  cyclones. 

Extensive  circumstantial  evidence  indicates  that  the  air- 
sea  fluxes  play  a  major  role  in  ocean  cyclogenesis .   Petterssen 
et  al. ,  (1962)  described  the  evolution  of  the  cyclone  based  on 
a  composite  of  51  cases  of  western  North  Atlantic  Ocean  storms. 
They  found  that  initial  storm  development  over  the  ocean  ap- 
pears to  be  due  mainly  to  low- level  thermal  advection.   This 
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is  in  contrast  to  development  over  land  where  upper-level 
vorticity  advection  is  the  dominant  development  mechanism. 
Nitta  and  Yamamoto  (1974)  studied  open  ocean  cyclogenesis 
over  the  East  China  Sea.   They  also  found  cases  of  cyclogene- 
sis occurring  with  zonal  flow  aloft,  which  implies  that  low- 
level  thermal  advection  is  the  dominant  development  mechanism 
over  the  ocean.   Other  case  studies  (e.g.,  Winston,  1955; 
Pyke ,  1965;  Sanders  and  Gyakum,  1980)  of  explosive  cyclogene- 
sis over  the  oceans  suggest  that  air-sea  exchanges  are  respon- 
sible for  the  observed  rapid  rates  of  development.   More 
indirect  evidence  has  been  obtained  recently  from  studies  of 
cyclogenesis  over  the  southern  oceans.   In  a  review  of  air- 
sea  interaction  in  the  southern  hemisphere,  Baker  (1979)  cites 
studies  that  show  the  highest  frequency  of  cyclone  development 
occurs  over  the  oceanic  polar  front.   Carleton  (19  81)  obtained 
the  same  result  using  satellite  data  to  locate  the  initial 
positions  of  developing  storms  in  the  southern  hemisphere. 
These  studies  are  only  able  to  postulate  that  air-sea  exchanges 
are  important  in  over-ocean  cyclogenesis.   Sanders  and  Gyakum 
(1980)  support  this  observation  and  also  find  a  positive  corre- 
lation between  the  intensity  of  the  cyclogenesis  and  the 
strength  of  the  ocean  temperature  gradient. 

Conflicting  evidence  regarding  the  importance  of  air-sea 
fluxes  has  been  reported  from  several  numerical  model  investi- 
gations and  one  observational  study.   Spar  and  Atlas  (1975) 
investigated  the  response  to  sea-surface  temperatures  using 
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the  Goddard  Institute  for  Space  Studies  general  circulation 
model.   They  found  the  model  to  be  insensitive  to  the  speci- 
fication of  sea-surface  temperatures  for  both  short-  and 
extended-range  forecasts.   Danard  and  Ellenton  (1979)  studied 
the  role  of  the  air-sea  fluxes  based  on  real  data  predictions 
with  an  eight-level  primitive  equation  model.   They  also  found 
that  the  air-sea  fluxes  are  unimportant  to  modelled  cyclo- 
genesis  along  the  east  coast  of  North  America.   Wei  (1979) 
presents  the  only  contrary  observational  study.   She  inves- 
tigated the  energetics  of  an  East  China  Sea  storm  using  the 
Air  Mass  Transformation  Experiment  (AMTEX)  data.   She  found 
that  surface  sensible  heat  fluxes  and  latent  heat  release  con- 
tribute very  little  to  the  total  available  potential  energy 
of  the  system,  and  concludes  that  the  diabatic  processes  are 
unimportant  in  ocean  cyclogenesis . 

The  mechanisms  by  which  air-sea  fluxes  could  enhance  ocean 
cyclogenesis  are  also  in  question.   Petterssen  et  al. ,  (1962) 
suggest  the  modification  of  low-level  thermal  advection  as  a 
mechanism,   pyke  (1965)  and  Simpson  (1969)  propose  an  indirect 
enhancement  by  increasing  cumulus  convection  which,  in  turn, 
increases  the  vorticity  of  the  cyclone.   Baker  (1979)  suggests 
that  ocean  temperature  gradients  induce  thermally  direct  mean 
meridional  circulations,  enhancing  the  jet  stream  in  their 
location  and  thereby  determining  the  region  of  cyclogenesis. 
Recent  studies  of  intermediate-scale  cyclones  developing  over 
the  ocean  reveal  similar  controversies.   Rasmussen  (1979) 
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analyzes  case  studies  of  polar  lows  over  the  eastern  North 
Atlantic  Ocean.   His  studies  indicate  that  these  lows  are 
convectively  driven,  implying  a  role  for  the  air-sea  fluxes 
similar  to  the  one  of  Pyke  (1965)  and  Simpson  (1969) .   Reed 
(1979)  and  Mullen  (1979)  presented  case  studies  of  similar 
polar  lows  in  the  eastern  North  Pacific  Ocean.   They  found 
the  lows  are  driven  by  baroclinic  instability  processes  with 
the  air-sea  fluxes  contributing  mainly  to  the  reduction  of 
low-level  static  stability  and  thereby  enhancing  development. 
Mullen  (1979)  and  Reed  (1979)  encourage  model  studies  to 
resolve  these  controversies. 

The  fundamental  hypothesis  of  this  study  is  that  the 
surface  fluxes  of  moisture  and  sensible  heat  significantly 
modify  the  planetary  boundary  layer  and  the  lower  atmosphere 
on  a  time  scale  of  one  to  four  days.   These  effects  change 
the  low-level  baroclinicity  and  static  stability  of  the  atmos- 
phere as  well  as  the  moisture  availability  for  latent  heat 
processes.   These  processes  may  account  for  the  observed 
differences  between  land  and  ocean  cyclogenesis . 

This  study  should  be  viewed  as  the  first  in  a  series  of 
numerical  experiments  at  the  Naval  Postgraduate  School  that 
are  designed  to  investigate  the  role  of  air-sea  interaction 
in  the  development  of  extratropical  cyclones.   Its  purpose 
is  to  examine  open-ocean  cyclogenesis  in  an  idealized  atmos- 
phere that  is  in  quasi-equilibrium  with  the  underlying  ocean- 
surface  temperatures.   The  study  does  not  include  a  land-sea 
boundary  and  is  not  intended  to  treat  the  more  energetic 
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cyclogenesis  which  occurs  near  the  coast,  presumably  as  a 
result  of  extreme  air-ocean  imbalances.   The  specific  objec- 
tives of  this  study  are  to: 

1.  determine  the  contribution  of  the  surface  fluxes  to 
the  growth  of  an  unstable  baroclinic  wave  during  the 
early  linear  growth  phase; 

2.  show  the  mechanisms  through  which  these  fluxes 
modify  the  growth  and  decay  of  a  baroclinic  wave; 
and 

3.  investigate  the  role  of  these  fluxes  in  the  develop- 
ment of  secondary  small-scale  cyclones. 

A  basic  assumption  of  this  research  is  that  the  scarcity 
of  conventional  observations  over  the  ocean,  and  the  complex 
interactions  of  the  physical  processes  involved,  require  the 
use  of  a  numerical  model.   In  particular,  the  model  should 
include  a  sophisticated  planetary  boundary  layer  parameteri- 
zation and  analytic  initial  conditions  to  avoid  uncertainties 
due  to  incomplete  or  incorrect  observations  over  the  oceans. 
Previous  observational  studies  (see  Pyke,  1965;  Laevastu, 
1965;  Simpson,  1969)  were  severely  hampered  by  a  lack  of  both 
conventional  observations  and  measurements  of  boundary  layer 
fluxes  required  to  investigate  these  interactions.   Previous 
numerical  studies  (Spar  and  Atlas,  1975;  Gambo,  1976;  Danard 
and  Ellenton,  1979)  used  either  conventional  observations  or 
numerical  models  in  which  the  parameterizations  of  the  dia- 
batic  processes  were  greatly  simplified.   This  experiment 
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utilizes  a  version  of  the  UCLA  general  circulation  model  that 
is  being  developed  by  the  Naval  Environmental  Prediction  Re- 
search Facility  (NEPP.F)  for  use  as  the  Navy  Operational  Global 
Atmospheric  Prediction  System  (NOGAPS) .   The  NOGAPS  model  has 
one  of  the  most  complete  parameterizations  of  the  interactions 
between  the  planetary  boundary  layer  and  the  free  atmosphere. 
This  model  provides  a  state-of-the-art  tool  for  investigating 
the  role  of  the  air-sea  fluxes.   However,  this  is  only  a 
model  study  and  must  be  validated  when  adequate  atmospheric 
observations  become  available. 

Chapter  II  of  this  paper  describes  the  NOGAPS  prediction 
model  used  in  this  research,  and  the  atmosphere  and  ocean 
conditions  chosen  to  initialize  the  model.   The  analysis  of 
the  experimental  results  is  divided  into  three  parts  for  ease 
of  interpretation.   These  are: 

1.  Chapter  III--the  adiabatic-diabatic  comparison  and 
the  large-scale  response  to  the  surface  fluxes, 

2.  Chapter  IV — the  initial  growth  phase  of  the  cyclone 
waves,  and 

3.  Chapter  V--the  development  of  the  secondary  cyclones. 
Chapter  VI  is  a  discussion  of  the  results  of  this  investiga- 
tion in  terms  of  the  three  objectives. 


19 


II.   THE  ATMOSPHERIC  MODEL 

The  model  chosen  for  this  study  is  a  state-of-the-art 
operational  forecast  model.   Extensive  evaluation  and  veri- 
fication of  the  model  has  been  performed  by  research  groups 
at  the  Naval  Oceanography  Command  (Paine,  1980)  and  at  UCLA 
(Mechoso  et  al_.  ,  1978)  .   The  model  has  been  chosen  by  the 
Naval  Oceanography  Command  to  replace  their  current  opera- 
tional atmospheric  prediction  model.   Any  model  is  an  ideali- 
zation of  the  real  world  limited  both  by  parameterizations 
of  real  world  processes  and  by  initial  conditions.   By  using 
a  numerical  model  in  this  study  it  is  hoped  that  the  results 
will  aid  in  the  interpretation  of  future  observational  studies 
and  also  provide  insight  into  the  difficulties  of  current 
forecast  models  in  handling  open-ocean  cyclogenesis  (Sanders 
and  Gyakum,  19  80) . 

The  following  sections  describe  the  NOGAPS  forecast  model, 
the  method  of  removing  diabatic  processes  from  the  model  and 
the  method  of  initializing  the  model. 

A.   FINITE  DIFFERENCE  SCHEME 

The  adiabatic  portion  of  the  model  is  described  in  detail 
by  Arakawa  and  Lamb  (1977),  and  thus  only  a  brief  review  is 
presented  here.   The  NOGAPS  code  for  this  study  was  kindly 
provided  by  Dr.  T.  Rosmond  of  the  Naval  Environmental  Pre- 
diction Research  Facility  (NEPRF)  and  contains  all  modifica- 
tions made  to  the  system  through  July,  1980.   A  smoothing  of 
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the  meridional  wind  component  on  the  latitude  rings  immedi- 
ately adjacent  to  the  poles  has  been  added  to  maintain  model 
stability  when  computing  on  the  fine  resolution  grid  (2°  lat 
by  2.5°  long) . 

The  prognostic  variables  for  the  model  are  the  horizontal 
velocity,  V,  temperature,  T,  a  form  of  the  surface  pressure, 
it  ,  and  specific  humidity,  q.   Additional  prognostic  variables 
associated  with  the  planetary  boundary  layer  (PBL)  will  be 
described  in  the  next  section.   The  model  uses  a  sigma  coor- 
dinate system  defined  as 


P  ~  Pj 


a      = (2-1) 


where: 

p    =   50  mb 
ir    =   Ps  -  Pr- 

p  is  pressure  and  p   is  the  pressure  at  the  earth's  surface. 
In  this  study  the  lower  boundary  of  the  model  is  open  ocean. 
The  top  of  the  model  atmosphere,  p  ,  is  50  mb.   No  strato- 
spheric physics  are  contained  in  this  version  (NOGAPS)  of  the 
model.   The  vertical  distribution  of  the  six  model  layers, 
assuming  an  arbitrary  surface  pressure  of  1000  mb ,  is  pro- 
trayed  in  Fig.  2-1.   The  vertical  velocity,  a,    is  carried  at 
the  layer  interfaces,  with  V,  T  and  q  defined  at  the  layer 
mid-levels.   Simmons  and  Hoskins  (1976)  investigated  the 
sensitivity  of  atmospheric  numerical  predictions  to  the  number 
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Figure  2-1.   Vertical  distribution  of  model  large-scale 
prognostic  variables.   Pressure  values  of 
sigma  levels  vary  with  surface  pressure.   A 
surface  value  of  1000  mb  is  assumed  in  this 
figure. 
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of  model  layers  in  the  vertical.   They  compared  results  from 
two-,  five-,  eight-  and  16-layer  versions  of  the  same  model. 
They  found  the  five-layer  model  accurately  represented  phase 
speeds  to  wave  number  16  (the  limit  of  their  study)  and  growth 
rates  to  wave  number  12.   For  waves  shorter  than  zonal  wave 
number  12,  the  five-layer  version  overestimates  growth  rates 
while  the  eight-layer  version  accurately  predicts  phase  speed 
and  growth  rates  to  the  limit  of  their  study,  wave  number  16. 
Six  layers  is  adequate  for  this  study,  however  a  greater  reso- 
lution in  the  lower  atmosphere  would  be  preferable  for  future 
studies  dealing  with  small-scale  waves. 

The  variables  are  staggered  in  the  horizontal  (Fig.  2-2) 
such  that  the  meridional  wind  component,  v,  is  carried  at 
points  north  and  south  of  the  center  point  and  the  zonal  wind 
component,  u,  is  carried  at  points  east  and  west  of  the  center 
point.   The  u  and  v  components  are  averaged  to  the  center  (T) 
point  for  diabatic  and  frictional  computations.   A  coarse- 
model  resolution  of  4°  lat  by  5°  long  was  used  during  initial 
model  development  to  reduce  integration  time.   This  resolution 
was  increased  to  2°  lat  by  2.5°  long  for  the  experiments  des- 
cribed below  to  properly  handle  growth  rates  and  phase  speeds 
of  the  large-scale  cyclones.   The  finer  resolution  results  in 
grid  distances  of  222  km  by  196  km  at  45°N/S,   which  is  also 
sufficient  to  resolve  secondary  systems  with  a  wavelength  of 
1500  to  2000  km. 

The  experiments  are  performed  on  a  60  degree  pole-to-pole 
sector  having  24  grid  points  east-west  and  91  grid  points 
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Figure  2-2. 


Horizontal  distribution  of  model  large- 
scale  prognostic  variables.   ua  which  is 
not  shown  is  carried  at  T-points . 
(Arakawa  scheme  C) 
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north-south.   The  model  uses  cyclic  east-west  boundary  con- 
ditions, which  allow  a  continuous  over-ocean  propagation. 
The  60  degree  sector  was  chosen  to  reduce  computation  time 
and  space  requirements.   Adiabatic  model  studies  by  Simmons 
and  Hoskins  (1978)  have  shown  simultaneous  development  of 
identical  waves  around  the  globe  when  the  model  was  initial- 
ized by  a  normal-mode,  wave  number  six  perturbation  as  in 
these  experiments.   Although  this  may  not  be  strictly  true 
for  this  diabatic  model  study,  the  differences  are  only  due 
to  small  time  differences  in  solar  heating  and  should  not 
affect  the  qualitative  aspects  of  the  model  experiments. 

The  model  utilizes  a  leap-frog  time  differencing  scheme 
with  a  time  step  of  three  minutes  in  the  fine  mesh  version. 
The  model  diabatic  package  is  executed  every  30  minutes,  and 
is  followed  by  a  Matsuno  time  step  to  control  the  computa- 
tional mode  and  to  assist  in  assimilation  of  the  diabatic 
effects  (see  Haltiner  and  Williams,  1980).   To  avoid  using 
an  extremely  short  time  step  to  maintain  computational  sta- 
bility near  the  poles,  a  smoothing  technique  is  used  on  the 
zonal  mass  flux  and  east-west  pressure  gradient  terms  in  that 
region  (Arakawa  and  Lamb,  1977).   An  additional  smoothing  of 
the  meridional  component  immediately  adjacent  to  the  poles 
is  required  for  the  fine-mesh  version. 

B.   DIABATIC  PARAMETERIZATION 

One  of  the  assumptions  of  this  research  is  that  a  sophis- 
ticated PBL  parameterization  (Table  2-1)  is  required  to 
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Table  2-1 
Summary  of  the  Model  Features 

Model  Dynamics 
60°  pole-to-pole  sector 
Ax   =   2.5°  long  (213  km  at  40°N/S) 
Ay   =   2.0°  lat  (240  km) 
Cyclic  east-west  boundary  conditions 
Staggered  grid  (Arakawa  scheme  C) 
a     vertical  coordinate 
6  model  layers  in  troposphere 

Model  Diabatics 
Planetary  boundary  layer  (Deardorf f/Randall) 

diagnostic  determination  of  surface  fluxes 

variable  depth 

stratus  cloud  layer  parameterization 
Arakawa-Schubert  cumulus  parameterization  (Lord) 
Large-scale  precipitation 
Moist  convective  adjustment 
Radiation  (Katayama/Schlesinger) 

diurnal  and  seasonal  variation 

cloud  and  water  modifications 
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adequately  simulate  the  direct  and  indirect  effects  of  the 
surface  fluxes  on  extratropical  cyclone  development.   The 
PBL  treatment  in  this  model  follows  Deardorff  (1972),  as 
implemented  in  the  UCLA  model  by  Randall  (1976) .   It  allows 
for  interactions  between  the  PBL  and  cumulus  cloud  ensembles 
and/or  a  stratus  cloud  layer  at  each  model  gridpoint.   Sur- 
face fluxes  of  momentum,  heat  and  moisture  are  diagnostically 
determined  using  a  bulk  Richardson  number  based  on  the  sea- 
surface  temperature  and  the  values  of  V,  T  and  q  from  the 
dynamic  portion  of  the  model.   These  values  are  then  used  to 
predict  a  new  PBL  depth  and  the  strength  of  the  inversion 
jumps  of  V,  T  and  q  at  the  top  of  the  PBL.   The  NOGAPS  PBL 
differs  from  that  described  in  Arakawa  and  Lamb  (1977)  in 
that  it  is  constrained  to  remain  within  the  lowest  model  layer, 
i.e.,  beneath  approximately  8  00  mb,  to  avoid  the  unrealistic 
PBL  depths  which  occurred  in  previous  versions  tested  at  NEPRF 
and  UCLA. 

Cumulus  parameterization  follows  the  scheme  of  Arakawa 
and  Schubert  (1974)  as  introduced  into  the  model  by  Lord 
(1978) .   The  base  of  the  cumulus  cloud  ensembles  is  taken  as 
the  top  of  the  PBL  and  provides  communication  between  the  PBL 
and  the  higher  model  layers  through  detrainment  of  mass  and 
moisture  at  the  cloud  top  level.   Condensation  also  occurs 
at  grid  points  where  the  air  becomes  supersaturated,  and  in 
a  moist  convective  adjustment  procedure  that  has  been  added 
to  remove  moist  convective  instability  between  layers  not 
associated  with  the  PBL.   The  radiational  heating  computation 
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which  follows  Katayama  (1972)  includes  the  effects  of  water 
vapor  and  also  cloud  distributions  predicted  by  the  PBL  and 
cumulus  parameterizations  and  by  large-scale  precipitation 
processes . 

C.   MODEL  VARIATIONS 

The  approach  used  in  these  experiments  is  to  selectively 
withdraw  the  physical  processes  from  the  complete  diabatic 
model.   This  results  in  much  less  model  shock  during  the  first 
few  hours  of  integration  than  in  experiments  where  physical 
processes  are  added  to  the  adiabatic  model  integration.   In 
the  adiabatic  case,  the  air-sea  differences  and  surface  winds 
can  become  abnormally  large.   This  can  result  in  excessive 
surface  heat,  moisture  or  momentum  fluxes  as  those  processes 
are  added  to  the  model.   In  the  diabatic  model,  the  air-sea 
fluxes  are  in  quasi-equilibrium  at  the  time  that  the  physical 
processes  are  withdrawn  from  the  model.   Thus,  there  is  al- 
most no  shock  generated  and  the  wave  cyclone  is  only  gradually 
modified. 

The  experiments  are  based  on  six  different  model  configura- 
tions.  The  nomenclature  to  be  used  in  describing  these  basic 
configurations  is 

1.  D — the  full  diabatic  model  as  described  above; 

2.  A — an  adiabatic  model  created  by  removing  the  physi- 
cal parameterization  package  which  contains  all 
diabatic  effects  and  friction,  and  by  adding  dry 
convective  adjustment  to  maintain  stability; 
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3  .   F — an  adiabatic  model  with  friction  represented  in 
the  PBL  parameterization,  but  not  allowing  any  heat 
or  moisture  transfer  to  the  large-scale  variables; 

4.  D-FS — the  full  diabatic  model  without  surface  sensi- 
ble heat  fluxes,  created  by  setting  the  surface  heat 
transfer  to  zero  within  the  PBL  calculations; 

5.  D-FQ — the  full  diabatic  model  without  surface  mois- 
ture fluxes,  created  as  in  (4)  but  for  the  moisture 
transfer; 

6.  D-FSQ — the  full  diabatic  model  without  surface  fluxes 
of  moisture  or  sensible  heat,  created  as  in  (4)  and 
(5)  . 

It  should  be  emphasized  that  this  is  intended  to  be  a  study 
of  physical  processes,  not  of  a  model  parameterization  scheme 
or  of  an  initialization  test  of  questionable  or  erroneous 
data.   These  experiments  start  with  a  system  in  which  the 
surface  fluxes  are  relatively  small.   Removing  the  surface 
heat  flux  does  not  remove  radiative  sources  and  sinks  of  heat- 
ing, nor  does  removing  the  surface  moisture  flux  prevent  con- 
densation of  or  the  redistribution  of  the  existing  moisture 
by  the  various  cloud  parameterizations  within  the  model. 
These  experiments  are  similar  to  "spin-down"  experiments 
and  as  such  they  require  time  before  the  absence  of  the  sur- 
face transfer  processes  is  felt. 

D.   INITIALIZATION 

The  model  initial  conditions  have  been  chosen  to  study 
the  air-ocean  fluxes  during  the  spring  and  fall  regimes  over 
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the  open  ocean.   The  model  start  time  has  been  selected  as  1 
October.   This  makes  the  northern  sector  a  fall  regime  and 
the  southern  sector  a  spring  (April)  regime  with  respect  to 
the  solar  fluxes . 

The  model  atmosphere  initial  conditions  are  analytically 
specified  following  Simmons  and  Hoskins  (1977,1978).   Zonal 
jets  are  centered  at  35°S  and  45°N  with  a  maximum  of  52  m/s 
in  the  southern  or  spring  hemisphere  and  3  5  m/s  in  the  northern 
or  fall  hemisphere  (see  Fig.  2-3) .   The  wind  profiles  follow 
the  function  u(a,cj>)  =  UQ(a)M(y)  where: 

<j)   =   latitude 

a   =   vertical  coordinate 

y   =   sin  <p 


M(y) 


sin3(7Ty2)       45°N  jet 


4    1247 
sin4  Uy   Z4/)   35°S  jet 


and  U.(a)  is  a  polynomial  function  approximation  of  the 
January  mean  zonal  wind  profile  at  30°N  taken  from  Oort  and 
Rasmussen  (1971).   This  polynomial  function  is  U0(cr)  for  the 
fall  profile  but  is  multiplied  by  a  factor  of  1.5  at  all 
levels  to  create  UQ(a)  for  the  spring  profile. 

The  jet  locations  and  intensities  are  intended  only  to 
represent  conditions  which  could  exist  during  those  seasons. 
April  jet  speeds  frequently  exceed  8  0  m/s,  the  maximum  at- 
tained in  these  experiments,  and  jet  latitudes  are  highly 
variable  (Palmen  and  Newton,  19  69) .   The  mean  latitudes  are 
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representative  of  spring  and  fall  conditions  as  given  in  Oort 
and  Rasmussen  (1971).   In  addition,  these  jets  guarantee  the 
growth  of  baroclinically  unstable  waves,  as  found  by  Simmons 
and  Hoskins  (1977,1978). 

Temperature  fields  are  derived  using  the  nonlinear  balance 
routine  developed  by  Mr.  E.  Barker  of  NEPRF.   The  mean  Janu- 
ary geopotential  heights  at  30 °N  from  Oort  and  Rasmussen 
(1971)  are  used  as  the  reference  profile  to  solve  for  the 
geopotential  field  from  the  above  analytic  wind  field  using 
the  nonlinear  balance  equation.   The  initial  temperature  field 
is  then  determined  from  the  geopotential  field  using  the  hydo- 
static  equation.   The  reference  geopotential  profile  is  speci- 
fied at  30°S  (spring  regime)  in  the  balance  routine  to  match 
the  temperature  field  as  closely  as  possible  to  the  clima- 
tological  sea-surface  temperature  gradient. 

In  order  to  avoid  large  initial  surface  fluxes,  emphasis 
>   has  been  placed  on  the  specification  of  low-level  temperature 
and  moisture  profiles.   In  these  experiments,  the  ocean  sur- 
face temperatures  are  independent  of  longitude  and  fixed  in 
time.   Ocean  surface  temperatures  are  derived  from  the  model 
surface  air  temperatures,  so  that  initially  the  lowest  model 
layer  is  near  neutral  stability.   Some  adjustments  are  made 
to  achieve  both  a  realistic  sea-surface  temperature  gradient 
and  a  field  that  is  consistent  with  the  simple  atmospheric 
structure  represented  by  the  specified  jet  intensity/location 
in  each  hemisphere.   In  Fig.  2-4  the  model  sea-surface  tempera- 
tures are  compared  with  smoothed  values  at  approximately  160 °W 
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taken  from  Robinson  (197  6) .   Climatological  values  for  the 
southern  hemisphere  (spring)  are  April  northern  hemisphere 
values.   The  lower  equatorial  temperatures  are  desirable  to 
reduce  the  influence  of  the  convectively-driven  Hadley  circu- 
lation.  The  specification  of  open-ocean  conditions  (no  sea 
ice)  is  responsible  for  the  temperature  differences  near  the 
poles. 

Realistic  moisture  profiles  for  the  marine  atmosphere  are 
more  difficult  to  specify.   The  intent  in  this  experiment  is 
to  represent  open-ocean  conditions,  i.e.,  mid-Pacific  or 
eastern  Atlantic.   Therefore,  it  is  natural  to  assume  much 
higher  relative  humidity  values  in  the  lower  model  levels  than 
the  global  mean  values.   Initial  moisture  fields  for  the 
model  were  specified  with  a  90%  relative  humidity  at  900  mb , 
and  a  linear  decrease  to  40%  at  the  highest  model  level  (70 
mb) .   Fig.  2-5  is  a  comparison  of  model  relative  and  specific 
humidities  with  zonal  mean  values  from  the  global  climatology 
by  Oort  and  Rasmussen  (19  71)  for  the  three  lowest  model  levels 
Model  relative  humidities  compare  favorably  with  climatology 
when  considered  as  open-ocean  versus  global  means.   Mullen 
(1979)  found  low-level  relative  humidities  of  80%  to  85% 
over  the  eastern  North  Pacific  Ocean.   Specific  humidities 
are  significantly  different  near  the  equator  and  the  poles 
due  to  differences  between  mean  atmosphere  and  model  tempera- 
tures in  those  regions. 

The  low-level  winds  are  initially  less  than  3  m/s .   This 
helps  to  prevent  an  initial  shock  or  imbalance  due  to  large 
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air-sea  fluxes.   Consequently,  the  PBL  is  set  to  a  very- 
shallow  depth  (10  mb)  with  no  temperature,  moisture  or 
velocity  jumps  at  the  top.   Rather  than  attempt  to  specify 
these  values,  it  is  intended  that  both  the  low-level  wind 
field  and  the  PBL  will  evolve  to  realistic  values  as  the 
initial  perturbation  develops. 

The  conditions  discussed  above  are  used  for  all  experi- 
ments with  unperturbed  initial  conditions  (hereafter  referred 
to  as  Z  for  zonal) ,  and  as  the  basis  for  creating  two  differ- 
ent perturbed  initial  conditions.   The  first  perturbed  case 
(S)  is  created  following  Simmons  and  Hoskins  (1977)  .   A 
small  amplitude  sinusoidal  perturbation  is  added  to  the 
meridional  wind  component  at  all  levels  prior  to  balancing. 
This  perturbation  is  in  zonal  wave  number  six  with  maximum 
values  of  .5  m/s  at  45°N/S,  and  equal  to  zero  at  the  equator 
and  the  poles.   The  associated  surface  pressure  perturbation 
is  approximately  1  mb.   These  initial  conditions  are  used  for 
the  extratropical  cyclone  evolution  experiments .   The  second 
perturbed  case  (R)  is  created  by  inserting  similar  small 
amplitude  perturbations  into  the  zonal  initial  conditions 
(Z)  after  the  fields  are  balanced.   In  this  case,  perturba- 
tions with  a  maximum  value  of  .1  m/s  in  wave  numbers  6,  12, 
18,  24,  30  and  3  6  are  added  to  the  meridional  wind  component. 
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III.   LARGE  SCALE  RESPONSE 

The  experimental  results  have  been  divided  into  three 
parts.   This  chapter  describes  the  large-scale,  baroclinic 
wave  growth  in  the  diabatic  model  atmosphere,  and  how  this 
compares  to  the  adiabatic  model  experiment.   The  adiabatic 
experiment  is  also  verified  with  a  similar  model  experiment 
by  Simmons  and  Hoskins  (1978).   This  chapter  introduces  the 
surface  flux  experiments  and  describes  their  results  in  terms 
of  large-scale  cyclone  development.   The  next  chapter  con- 
siders the  effects  of  the  individual  diabatic  processes  dur- 
ing the  linear  growth  phase  (first  72  h)  of  the  baroclinically 
unstable  waves.   In  the  following  chapter,  the  development 
of  secondary  cyclone  waves  in  the  diabatic  model  experiment 
is  discussed.   In  particular,  the  effect  of  the  surface  fluxes 
of  sensible  heat  and  moisture  on  the  development  of  the  sec- 
ondary cyclone  waves  is  examined  based  on  these  model 
experiments . 

The  experiments  are  of  two  basic  types--those  starting  at 
day  0,  and  those  initialized  at  day  4  using  fields  from  the 
full  diabatic  model  integration  (Table  3-1) .   The  objective 
of  the  first  group  is  to  establish  the  magnitude  of  the  model 
atmosphere  response  to  the  diabatic  parameterizations .   The 
second  group  of  experiments  is  performed  to  study  the  influ- 
ence of  the  individual  diabatic  processes.   These  experiments 
employ  the  technique  (described  in  Section  II. C)  of  selectively 
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Table  3-1 
Experiments  for  the  Large-scale  Study 

*  * 

Exp.     Initial  Integration      Model 

Fields  Period         Version 

Experiments  Initialized  from  Day  0 
D        S  (wave  no.  6)       15  days       complete  diabatic 
AS  15  days       adiabatic 

DZ       Z-zonal  3  days       complete  diabatic 

DR       R-multiwave  5  days       complete  diabatic 

Experiments  Initialized  from  Day  4 

F        S4  (fields  from      6  days       adiabatic  with 
D  at  day  4)  model  friction 

D-FS     S4  6  days       diabatic  w/o  surface 

sensible  heat  flux 

D-FQ      S4  6  days        diabatic  w/o  surface 

moisture  flux 

D-FSQ    S4  6  days       diabatic  w/o  surface 

moisture  or  sensible 
heat  fluxes 


* 


See  Sections  II. C  and  II. D  for  a  more  complete  discussion 
of  model  variations  and  initial  conditions. 
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removing  the  model  parameterizations  to  investigate  the 
contributions  of  the  surface  fluxes  of  sensible  heat  and 
moisture. 

A.   DIABATIC-ADIABATIC  COMPARISON 

The  adiabatic  model  experiment  (A)  compares  very  closely 
with  the  results  of  a  similar  study  by  Simmons  and  Hoskins 
(1978)  .   Surface  pressure  maps  for  experiment  A  at  day  6  are 
displayed  in  Figs.  3-1  a-b.   The  northern  hemisphere  case 
(3-la) ,  which  employs  initial  conditions  adapted  from  Simmons 
and  Hoskins  (1978),  is  essentially  the  same  as  in  their  study 
(see  their  Fig.  1) .   The  surface  pressure  field  exhibits  a 
pure  wave  number  six  pattern  in  each  hemisphere.   The  growth 
of  the  cyclone-anticyclone  pair  is  much  greater  in  the  spring 
(southern)  hemisphere,  as  would  be  expected  with  the  more  in- 
tense initial  jet  in  that  hemisphere.   Southern  hemisphere 
maps  in  Fig.  3-lb  and  below  are  presented  in  inverted  form 
(south  at  top,  east  remains  toward  right)  for  ease  of  comparison 
with  the  northern  hemisphere  results.   The  positions  of  the 
cyclone  centers  also  reflect  the  more  equatorward  initial 
latitude  of  the  upper-level  jet  maximum  in  the  spring  hemi- 
sphere.  The  vertical  structure  of  the  growing  cyclone  (Fig. 
3-3)  is  that  of  a  typical  growing  baroclinic  wave  as  seen  in 
observational  studies  (Palmen  and  Newton,  1969) .   This  struc- 
ture also  follows  closely  that  obtained  by  Hoskins  (1978) 
in  a  theoretical  study  of  baroclinic  waves.   In  agreement 
with  Simmons  and  Hoskins  (1978),  the  wave  cyclone  in  A  reaches 
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Figure    3-1. 


Surface  pressure  for  A  after  six  days  of 
integration  for   a.  northern  hemisphere  and 
b.  southern  hemisphere.   Contour  interval 
4  mb .   Southern  hemisphere  inverted  with  pole 
(south)  at  top  (east  remains  toward  right) . 


30 
LONGITUDE 


Figure  3-2.   As  in  Fig.  3-1  except  for  diabatic  model 
(D)  . 
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a  maximum  intensity  at  about  day  9  and  then  decays  slowly 
through  day  15  of  the  model  experiment.   Simmons  and  Hoskins 
(1978)  demonstrated  that  this  decay  was  brought  about  by 
barotropic  processes  which  became  dominant  during  the  later 
part  of  the  forecast  period. 

The  full  diabatic  model  experiment,  D,  is  used  for  adia- 
batic-diabatic  comparisons  and  as  the  control  case  for  the 
surface  flux  model  experiments.   The  inclusion  of  diabatic 
processes  in  this  experiment  results  in  much  more  rapid  de- 
velopment of  the  primary  low  by  day  6  (see  Figs.  3-2  a-b 
and  3-1  a-b).   In  addition,  shorter  wavelength  disturbances 
grow  in  the  diabatic  experiment  as  found  by  Haltiner  (1967). 
Although  shorter  wavelengths  are  not  apparent  at  day  6  in 
the  spring  case  (Fig.  3-2b) ,  they  are  present  at  earlier 
(see  Fig.  3-6b  below)  and  later  times  (see  Chapter  V)  in 
the  diabatic  model  run.   Secondary  cyclone  waves  are  more 
prevalent  in  the  fall  hemisphere,  but  they  are  an  important 
feature  of  both  hemispheres  in  the  diabatic  model  experiments 

Experiments  A  and  D  are  initialized  with  a  wave  number 
six  perturbation.   Following  Simmons  and  Hoskins  (1978)  and 
other  studies  of  baroclinic  wave  growth,  it  was  assumed  that 
wave  number  six  would  be  the  most  unstable  for  these  experi- 
ments.  Instead,  the  initial  period  of  one  to  three  days  is 
dominated  by  wave  number  18  growth.   As  mentioned  above,  a 
wave  number  6  disturbance  eventually  evolves.   To  estimate 
the  sensitivity  to  the  use  of  wave  number  6  in  the  initial 
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conditions,  a  second  diabatic  experiment  (DR)  was  performed 
using  the  multi-wave  initial  conditions  described  in  Section 
II. D.   Wave  number  18  again  dominates  the  early  growth  period, 
but  a  distinct  wave  number  6  pattern  evolves  by  day  5  (Figs. 
3-4  a-b) .   Haltiner  (1967)  conjectured  that  the  evolution  of 
large-scale  waves  in  the  atmosphere  was  more  the  result  of  a 
change  in  scale  of  shorter  wavelengths,  than  the  growth  from 
a  small  amplitude  perturbation  of  large-scale  dimensions. 
This  is  apparently  true  in  D  and  DR.   This  aspect  of  the  dia- 
batic model  experiments  and  the  growth  of  the  smaller  scale 
waves  is  discussed  in  detail  in  the  following  chapters. 

The  primary  trough  in  D  exhibits  a  larger  phase  speed 
than  that  of  the  primary  trough  in  A  during  the  period  of 
strongest  development  (Figs.  3-5  a-b).   This  speed  differ- 
ence is  approximately  5°  long  per  day  (450  km  at  35°S)  at 
its  maximum;  however,  the  phase  speeds  are  approximately 
the  same  after  day  5. 

The  release  of  latent  heat  narrows  and  strengthens  by 
several  fold  the  zone  of  upward  motion.   This  increases  the 
conversions  of  eddy  available  potential  energy  to  eddy  kinetic 
energy,  and  enhances  the  growth  of  the  baroclinic  wave  (Palmen 
and  Newton,  1969;  Haltiner,  1971) .   This  more  rapid  growth 
is  evident  in  D.   A  result  not  anticipated  from  previous 
work  is  that  the  cyclones  in  A  eventually  (about  day  9)  be- 
come nearly  as  intense,  98  0  mb/980  mb  (northern/southern 
hemisphere)  central  pressure  versus  982  mb/976  mb,  as  in  D. 
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Figure  3-4 . 


Surface  pressure  for  DR  after  five  days 
of  integration  for   a.  northern  hemisphere 
and   b.  southern  hemisphere.   Contour 
interval  4  mb. 
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Figure  3-5 


Position  of  streamline  troughs  at  sigma 
level  3  (300  mb)  in  the  southern  hemisphere 
for   a.  A  and   b.  D. 
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Figure  3-6.   As  in  Fig.  3-2  except  after  four  days  of 
integration. 
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Beyond  day  8,  the  primary  lows  begin  a  gradual  decay  in 
intensity  as  in  A.   The  barotropic  decay  process  appears  to 
control  the  diabatic  experiment  as  in  the  adiabatic  experi- 
ment.  Thus,  although  the  minimum  surface  pressure  appears 
to  depend  on  the  large-scale  baroclinic  and  barotropic  proc- 
esses, the  diabatic  processes  greatly  affect  the  timing  and 
rate  of  development  of  the  primary  lows. 

in  summary,  the  adiabatic  model  experiment  is  consistent 
with  previous  studies,  notably  Simmons  and  Hoskins  (1978). 
The  primary  cyclones  attain  approximately  the  same  maximum 
intensity  in  the  diabatic  as  in  the  adiabatic  experiments 
and  undergo  a  similar  decay  period.   The  growth  of  the  pri- 
mary cyclone  is  much  more  rapid  in  the  diabatic  experiment, 
however,  and  involves  interactions  with  wave  scales  not  present 
in  the  adiabatic  experiment.   Also,  in  the  adiabatic  experi- 
ment the  hemispheric  differences  are  primarily  in  the  inten- 
sity of  the  primary  cyclone.   The  diabatic  experiment  demon- 
states  a  stronger  tendency  toward  shorter  wavelengths  in  the 
fall  hemisphere  as  well.   These  results  agree  closely  with 
cbservational  studies.   Sanders  and  Gyakum  (1980)  found  a 
greater  frequency  of  rapidly  deepening  storms  over  the  oceans 
where  diabatic  processes  are  more  dominant.   Also,  there  is 
a  greater  tendency  for  shorter  wavelength  cyclones  to  develop 
over  the  mid-  and  high-latitude  oceans  than  over  land  (Reed, 
1979;  Nitta  and  Yamamoto,  19  74) . 
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B.   SURFACE  FLUX  MODIFICATION  TO  THE  PRIMARY  LOW  DEVELOPMENT 
Four  model  versions  are  initialized  with  fields  from  day 
4  of  D  (Figs.  3-6  a-b)  to  study  the  contributions  of  the 
diabatic  processes  (Table  3-1) .   The  broad  areas  of  low  pres- 
sure on  the  right  side  of  each  figure  evolved  into  intense 
cyclones  (see  Figs.  3-2  a-b)  during  the  next  two  days  of 
integration  in  all  model  versions.   Removal  of  the  surface 
fluxes  produced  relatively  minor  changes  in  the  intensity  of 
the  primary  cyclones  and  anticyclones  (Table  3-2)  and  very 
slight  changes  in  phase  speed  (not  shown)  by  day  6,  which  is 
two  days  after  removal  of  the  fluxes.   Removal  of  the  surface 
sensible  heat  flux  (D-FS)  leads  to  an  additional  1-3  mb  deep- 
ening of  the  primary  cyclone.   The  largest  effect  occurs  in 
the  northern  or  fall  hemisphere  with  the  warmer  SST  and  weaker 
upper-level  jet.   The  removal  of  the  surface  moisture  flux 
(D-FQ)  leads  to  4-5  mb  less  deepening  of  the  primary  cyclone 
in  two  days.   No  further  deepening  occurred  beyond  day  6  in 
D-FQ  and  D-FSQ.   Thus,  there  are  differences  of  9-10  mb  in 
central  pressure  by  day  8,  four  days  after  initialization. 
These  minor  changes  are  expected  from  previous  observational 
studies  (see  Palmen  and  Newton,  19  69)  and  also  from  numerical 
model  studies  (Haltiner,  1971;  Spar  and  Atlas,  1975). 

The  surface  pressure  fields  at  day  6  of  the  no  surface 
flux  experiments  (Figs.  3-7  a-c)  are  very  similar  to  those 
of  the  diabatic  experiment  (Fig.  3-2a) .   However,  large  sur- 
face fluxes  of  moisture  (Fig.  3-8)  which  are  present  in  D 
and  D-FS,  but  are  not  present  in  D-FQ  and  D-FSQ,  lead  to 
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Table  3-2 

Central  pressures  (mb)  for  the  primary  cyclone/ 
anticyclone  by  experiment  for  days  4,  6  and  8. 


Model 
Exp. 

Hemisphere 

Time 
Day  4 

Day  6 

Day  8 

D 

Fall 

1012/1020 

988/1028 

983/1024 

Spring 

1002/1020 

975/1022 

981/1024 

D-FS 

Fall 

initialized 
with  fields 
from  D  above 

987/1027 

980/1022 

Spring 

ii 

974/1024 

980/1028 

D-FQ 

Fall 

initialized 
with  fields 
from  D  above 

992/1024 

992/1024 

Spring 

n 

979/1022 

985/1024 

D-FSQ 

Fall 

initialized 
with  fields 
from  D  above 

992/1024 

993/1024 

Spring 

ii 

978/1024 

985/1024 

F 

Fall 

initialized 

997/1029 

992/1028 

with  fields 
from  D  above 

Spring        "  997/1026     981/1025 
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Figure  3-7. 


Northern  hemisphere  surface  pressure  at 
day  6  for   a.  D-FS,   b.  D-FQ  and   c.  D-FSQ 
Contour  interval  4  mb. 
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Figure  3-8. 


Northern  hemisphere  surface  moisture  flux 
rate  at  day  6  for  D.   Contour  interval 
4  mm/day. 
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different  cyclone  evolutions  in  these  experiments.   The 

largest  moisture  fluxes,  which  are  equivalent  to  approxi- 

2 

mately  1200  cal/cm  /day  of  latent  heat,  are  located  in  the 

southwest  quadrant  of  the  primary  cyclone  in  D  and  in  D-FS 
(similar  to  D  but  not  shown) .   These  values  are  consistent 
with  previous  observational  studies  (e.g.,  Petterssen  et  al. , 
1962  and  others) .   After  two  days,  the  loss  of  the  moisture 
flux  has  led  to  large  changes  in  the  distribution  of  low- 
level  moisture  (Figs.  3-9  a-d) .   Very  dry  air  is  advected 
into  the  primary  cyclone  in  D-FQ  and  D-FSQ  (Figs.  3-9  c-d) 
as  a  result  of  the  lack  of  the  surface  moisture  flux.   This 
results  in  a  decrease  in  latent  heat  release  in  the  primary 
cyclone  and  a  large  reduction  in  its  growth  after  day  6 
(Table  3-2) . 

It  is  generally  accepted  that  the  temperature  gradient 
associated  with  the  cold  front  of  a  developing  cyclone  is 
greater  than  that  of  the  warm  front.   For  the  low-level 
temperature  fields  of  the  experiments  performed  here  (Figs. 
3-10  a-e)  this  is  only  true  in  the  adiabatic  experiment,  F, 
and  possibly  the  no  surface  flux  experiment,  D-FSQ.   Both  the 
surface  sensible  heat  flux  and  the  surface  moisture  flux  re- 
duce the  strengthening  of  the  low-level  temperature  gradient 
by  thermal  advection.   D-FS  and  D-FQ  have  a  larger  tempera- 
ture gradient  in  the  region  of  the  primary  cyclone  than  in  D. 
This  gradient  is  even  larger  in  D-FSQ,  and  approaches  that 
of  the  adiabatic  experiment,  F.   Although  this  does  not 
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Figure  3-9. 


Northern  hemisphere  specific  humidity  fields 
at  sigma  level  6  (900  mb)  at  day  6  for 
a.  D,   b.  D-FS,   c.  D-FQ  and   d.  D-FSQ. 
Contour  interval  2  g/kg . 
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Figure    3-10. 


As  in  Fig.  3-9  except  for  temperature 
fields.   Contour  interval  5°C.   Also 
shown  (next  page)   e.  experiment  F. 
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apparently  affect  the  evolution  of  the  primary  cyclone  on  a 
two-to-three  day  time  scale,  it  could  have  significant 
effects  on  the  development  of  secondary  or  subsequent 
cyclones . 

Diabatic  processes  play  a  greater  role  in  the  growth  of 
the  primary  cyclone  in  the  fall  hemisphere  with  the  weaker 
jet  and  higher  SST.   Over  the  open  ocean,  the  surface  fluxes 
of  moisture  and  sensible  heat  only  weakly  modify  the  develop- 
ment of  the  primary  cyclones  on  a  time  scale  of  less  than  two 
days,  as  found  in  the  case  studies  of  Mullen  (1979) .   The 
modifications  of  the  low-level  temperature  and  moisture  dis- 
tributions are  significant,  however.   In  particular,  the  sur- 
face fluxes  tend  to  weaken  the  temperature  and  moisture 
gradients  associated  with  strong  advection,  especially  in 
the  western  sector  of  the  cyclone  where  the  surface  fluxes 
are  largest  (Petterssen  et  al .  ,  1962) .   This  results  in  more 
intense  gradients  in  the  eastern  sector  than  in  the  western 
sector  of  the  primary  cyclone  in  the  model  marine  atmosphere. 
This  is  in  contrast  to  that  found  in  the  adiabatic  experi- 
ments or  in  typical  over-land  examples  of  cyclogenesis .   Be- 
yond two  days,  these  modifications  affect  the  development  of 
the  primary  cyclone.   The  response  of  secondary  or  subsequent 
cyclones  may  also  be  affected,  as  presented  in  Chapter  V. 
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IV.   INITIAL  DEVELOPMENT 

A  major  feature  of  the  diabatic  experiment  is  the  evolu- 
tion of  multiple  lows  and  secondary  lows  (see  Fig.  3-2a) .   A 
cyclone  scale  initially  much  shorter  than  wave  number  six 
expected  from  the  adiabatic  experiment  predominates  during 
the  first  three  to  four  days  of  the  diabatic  experiment. 
Beyond  day  4,  as  seen  in  Figs.  3-6  and  3-2,  these  multiple 
lows  gradually  evolve  into  an  intense  primary  low  with  asso- 
ciated secondary  wave  cyclones.   The  evolution  of  these  second- 
ary cyclones  is  addressed  in  Chapter  V.   This  chapter  describes 
the  initially  linear  growth  phase  and  the  reasons  why  the 
smaller-scale  cyclones  occur. 

Table  4-1  lists  the  model  experiments  viewed  in  this 
chapter.   The  initial  conditions  and  model  configurations 
used  for  these  experiments  are  presented  in  Chapter  II. 

Table  4-1 

Experiments  for  the  initial  development  study 

Exp.     Initial  Integration         Model 

Fields  Period  Version 

D       S — wave  no.  6      15  days       complete  diabatic 

AS  15  days       adiabatic 

A36      Si. 5  (fields  from    2  days        adiabatic 
D  at  day  1.5) 

D-FS     S  3  days        diabatic  w/o  surface 

sensible  heat  flux 

D-FQ     S  3  days        diabatic  w/o  surface 

moisture  flux 

* 
See  sections  II. C  and  II. D  for  a  more  complete  discussion 

of  model  variations  and  initial  conditions. 
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A.   DIABATIC  CONTRIBUTIONS  TO  WAVE  GROWTH 

Growth  rates  by  wavenumber  at  300  and  900  mb  for  experi- 
ments A,  A36  and  D  are  displayed  in  Figs.  4-1  a-d.   The 
growth  rates  are  obtained  by  performing  a  Fourier  analysis 
of  the  meridional  wind  components  for  two  time  periods.   The 
exponential  growth  rate  for  each  wave  number  is  then  obtained 
from 


-1     V£k(t2} 
e   =   (t9-t,)  x  ln(    ,.A)  ,  (4-1) 

z   x       v£k^l; 


where: 

e  =  growth  rate 

v  =  meridional  wind  component 

k  =  wavenumber  (6,12,18,24) 

l  =  model  sigma  level ( -  900  mb  and  300  mb) 

t,  =  initial  time,  and 

ty  =  final  time. 

The  most  obvious  difference  between  experiments  A  and  D 
is  that  only  wave  number  6  grows  in  A.   This  has  already  been 
described  and  is  most  clearly  seen  in  the  surface  pressure 
fields  displayed  in  Figs.  3-1  and  3-6.   Experiment  A36  was 
included  to  incorporate  any  changes  in  the  model  basic  flow 
that  may  have  occurred  due  to  slight  imbalances  in  the  initial 
conditions,  or  due  to  the  adjustment  of  the  diabatic  model  to 
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Figure  4-1 


Exponential  growth  rate  between  days  2.5  and 
3.5  for  the  meridional  velocity  perturbation 
by  wavenumber  (see  text  for  description)  for 
a.  sigma  level  3  (300  mb) ,  36°S/   b.  sigma 
level  2,    46°N,   c.  sigma  level  6  (900  mb) , 
36°S  and   d.  sigma  level  6,  46°N.   Experi- 
ments D,  A  and  A36  are  compared. 
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the  underlying  sea-surface  temperatures.   Experiments  A  and 
A36  use  the  same  adiabatic  version  of  the  model—the  only 
difference  is  that  A36  is  initialized  using  fields  from  day 
1.5  of  the  diabatic  experiment  while  A  uses  the  day  0  fields. 
In  A36,  wave  number  12  grows  more  rapidly  than  wave  number  6 
in  the  lower  atmosphere  (9  00  mb) .   At  higher  levels,  wave 
number  12  exhibits  no  growth  in  A36.   Another  readily  evident 
feature  is  the  north/south  variation  in  the  A3 6  experiment  at 
higher  wave  numbers.   Wave  numbers  18  and  24  grow  weakly  at 
low  levels  and  not  at  all  in  upper  levels  in  the  weaker  jet/ 
higher  SST  fall  hemisphere.   In  the  spring  hemisphere,  these 
wave  numbers  exhibit  a  growth  comparable  to  wave  number  6  at 
both  levels.   The  last  feature  to  note  relative  to  A36  is 
that  wave  number  6  grows  more  rapidly  in  this  experiment  than 
in  A. 

The  diabatic  experiment  D,  exhibits  a  definite  shift  of 
the  most  unstable  wave  to  shorter  wavelengths.   At  upper 
levels,  wave  number  12  is  the  most  rapidly  growing  wave  in  both 
hemispheres,  while  at  low  levels  wave  number  18  dominates. 
This  is  a  distinctly  different  result  from  that  of  the  adia- 
batic experiments.   Wave  number  6  also  grows  more  rapidly  in 
D  than  A  everywhere  except  at  low-levels  in  the  northern  or 
fall  hemisphere.   Figs.  4-2  a-b  give  the  streamline  troughs 
versus  height  and  longitude  for  A  and  D.   These  figures  pro- 
vide a  clearer  picture  of  the  actual  wave  number/height  rela- 
tionships described  in  the  growth  rate  calculations,  although 
variations  in  intensity  are  not  indicated. 
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Figure  4-2 


Location  of  streamline  troughs  at  day  3 
for   a.  A  and   b.  D. 
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Several  model  studies,  e.g.,  Staley  and  Gall  (1977)  and 
Blumen  (1979),  have  found  that  the  wavelength  of  the  most 
unstable  wave  in  the  atmosphere  is  very  sensitive  to  changes 
in  low-level  static  stability.   In  addition,  Rao  and  Ferreira 
(1979)  performed  a  numerical  study  of  intermediate-scale 
(less  than  2000  km)  waves  using  a  moist  model  with  a  CISK- 
type  parameterization  of  convective  heating.   They  found  that 

the  use  of  a  standard  atmosphere  value  for  static  stability 

2  -2   -2 
of  . 02  m  s   mb   produced  two  modes  of  unstable  waves.   The 

first  was  a  conventional,  long  baroclinic  wave  in  which  latent 
heating  was  relatively  unimportant.   The  second  was  an  inter- 
mediate-scale disturbance  in  which  latent  heating-induced 
vertical  velocities  were  of  the  same  order  or  larger  than 

those  of  the  baroclinic  processes.   When  a  small  value  of 

2  -2   -2 
static  stability  (.002  m  s   mb   )  was  used,  only  a  conventional 

baroclinic  wave  was  found,  but  with  an  intermediate-scale 

wavelength. 

The  above  study  and  a  similar  one  by  Gambo  (1976)  suggest 

that  there  are  two  modes  of  intermediate-scale  unstable  waves: 

1.  a  latent  heat-driven  mode  with  energetics  similar  to 
a  tropical  system  but  not  dependent  on  low  values 

of  static  stability;  and, 

2.  a  conventional  baroclinic  mode  which  arises  due  to 
very  low  values  of  static  stability. 

Clearly  only  the  conventional,  baroclinic  mode  can  occur  in 
the  adiabatic  experiments.   While  the  diabatic  model  does  not 
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have  a  CISK  latent  heat  parameterization,  the  more  complete 
parameterizations  of  convective,  PBL  and  large-scale  latent 
heating  processes  should  allow  the  formation  of  either  inter- 
mediate-scale mode.   In  the  discussion  below,  differences  in 
wave  growth  among  experiments  A,  A36  and  D  are  explained  in 
terms  of  these  two  unstable  modes. 
Static  stability  is  defined  as 


a   3_9_ 
9   3p 


(4-2) 


where: 

a   =   specific  volume 

9   =   potential  temperature  and 

p   =   pressure. 

Figs.  4-3  a-b  give  the  variation  of  3  9/»p  as  a  function  of 
time  for  experiments  A,  A36  and  D  for  the  900-500  mb  layer. 

The  corresponding  values  of  static  stability  range  from  .015 

2  -2   -2 
to  .012  m  s   mb   .   Values  of  static  stability  are  computed 

separately  for  the  900  mb-surface  (D  only) ,  and  900-700  mb  and 

700-500  mb  layers.   The  lowest  values  are  found  in  the  900- 

2  -2   -2 
700  mb  layer  (.011  to  .013  m  s   mb   ) ,  but  the  greatest  tem- 
poral variability  is  found  in  the  700-500  mb  layer  (.015  to 

2  -2   -2 
.022  m  s   mb   ) .   The  averaging  of  the  two  model  layers 

seems  to  describe  the  variability  of  static  stability  in 

the  lower  model  atmosphere  better. 
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Figure    4-3. 


Change  in  potential  temperature  differences 
with  pressure  between  sigma  levels  4  and  6 
(500  to  900  mb)  in  °C/100  mb  for   a.  36°S 
and   b.  46°N. 
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In  the  region  of  the  jets  in  both  hemispheres,  the 
zonally-averaged  static  stability  reduces  rapidly  during  the 
first  two  to  three  days  of  integration  in  D.   The  trend  then 
reverses  around  day  3  to  3.5,  and  the  static  stability  in- 
creases to  a  value  comparable  to  A  by  day  5.   Both  experi- 
ments A  and  A36  maintain  a  relatively  constant  static  stability 
throughout  the  period  being  considered,  with  the  A3 6  case  hav- 
ing  a  much  smaller  value  due  to  the  initialization  at  day 
1.5  from  D.   It  is  important  to  note  that  the  value  of  static 
stability  is  similar  in  experiments  A36  and  D,  particularly 
in  the  southern  or  spring  sector. 

Low-level  static  stability  (Fig.  4-3)  and  the  growth  rates 
of  wave  numbers  6  and  12  (Fig.  4-1)  appear  to  be  associated, 
especially  in  the  adiabatic  experiments.   The  only  energy 
source  for  the  waves  in  the  adiabatic  experiment  is  baro- 
clinic  instability.   No  changes  in  static  stability  due  to 
air-sea  fluxes  or  latent  heat  release  can  occur  in  the  adia- 
batic experiments.   Thus,  it  appears  significant  that  wave 
number  12  develops  in  A36  at  low-levels,  but  does  not  develop 
in  A  which  has  a  larger  static  stability  (by  a  factor  of  40%) . 
Also,  Mullen  (1979)  found  that  small-scale  cyclones  in  the 
north  Pacific  Ocean  generally  have  very,  small  values  of  static 
stability.   As  noted  in  Gambo  (1976)  and  Rao  and  Ferreira 
(1979) ,  the  intermediate-scale  baroclinically  unstable  mode 
is  a  shallow  wave,  and  it  is  not  present  at  300  mb  in  A36. 
The  deeper  and  longer  (wave  number  6)  wave  is  relatively  in- 
sensitive to  low-level  static  stability  changes  and  its  growth 
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rate  in  A36  increases  only  slightly  over  that  of  wave  number 
6  in  A.   The  variations  in  growth  rates  of  wave  numbers  18 
and  24  at  the  900  mb  level  in  A36  may  be  similarly  explained 
in  terms  of  low-level  static  stability.   However,  the  growth 
of  these  waves  at  300  mb  (southern  case  only)  is  nearly  iden- 
tical to  that  of  wave  number  6  in  A36.   Because  these  waves 
(18  and  24)  are  expected  to  be  shallow  waves,  their  occurrence 
at  300  mb  is  likely  due  to  nonlinear  effects  (Williams,  1965) . 

Low-level  static  stability  in  D  is  not  significantly  dif- 
ference from  that  of  A36.   Nevertheless  the  wavelength  of 
the  fastest  growing  wave  is  shifted  to  shorter  wavelengths 
in  D  at  all  levels.   Wave  number  12  dominates  at  300  mb  and 
wave  number  18  dominates  at  90  0  mb.   However,  wave  number  6 
grows  only  slightly  faster  in  D  than  in  A.   Static  stability 
is  not  likely  to  be  the  cause  for  this  additional  shift,  and 
the  only  model  differences  are  in  the  diabatic  processes.   This 
implies  a  latent  heating  process  is  responsible  for  the  en- 
hanced growth  of  the  shorter  wavelength  waves,  in  agreement 
with  the  second  unstable  mode  of  Gambo  (19  76)  and  Rao  and 
Ferreira  (1979) .   Another  fact  supporting  this  reasoning  is 
the  much  larger  growth  rates  at  300  mb  in  the  northern  or  fall 
case.   In  terms  of  baroclinic  instability  processes,  the  much 
stronger  baroclinically  unstable  jet  is  located  in  the  southern 
case.   Consequently,  the  largest  growth  rates  should  be  pre- 
sent there,  which  is  true  for  the  shorter  waves  of  A36.   If, 
however,  latent  heat  release  is  responsible  for  the  larger 
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growth  rates,  the  lower  vertical  wind  shear  and  higher  sea- 
surface  temperatures  in  the  northern  case  should  produce  the 
largest  growth.   Also,  as  is  shown  below,  wave  number  18  com- 
mences rapid  growth  at  the  onset  of  D,  when  the  low-level 
static  stability  is  comparable  to  that  in  A.   This  growth 
appears  to  be  in  response  to  low-level  precipitation  or 
latent  heat  release  occurring  due  to  the  very  moist  initial 
atmosphere,  and  also  due  to  surface  heat  fluxes. 

Fig.  4-2b  also  depicts  a  more  vertical  structure  asso- 
ciated with  the  shallower  wave  in  D  than  in  the  deep  waves, 
or  in  wave  number  6  (Fig.  4-2a)  of  A.   The  east-west  tilt  of 
the  deeper  waves  is  more  typical  of  a  growing  baroclinic  wave, 
while  the  lack  of  tilt  of  wave  number  18  is  more  suggestive 
of  a  latent  heat-driven  wave  (Palmen  and  Newton,  1969) . 

To  summarize  this  section,  much  shorter  wavelength  waves 
are  present  in  the  diabatic  experiment  than  in  the  adiabatic 
experiments.   Two  processes  seem  to  be  responsible  for  the 
growth  of  the  shorter  wavelength  waves.   The  first  is  the 
reduction  of  the  low-level  static  stability  by  the  diabatic 
processes.   The  second  is  the  release  of  latent  heat.   The 
growth  rates  in  the  fall  case  with  higher  sea-surface  tempera- 
tures and  a  weaker  jet  responded  to  these  diabatic  effects 
more  strongly  than  in  the  spring  case  with  the  stronger  baro- 
clinically  unstable  jet.   These  results  are  in  close  agreement 
with  the  case  studies  of  Reed  (1979)  and  Rasmussen  (1979). 
Reed  observed  small  values  of  low-level  static  stability 
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associated  with  the  development  of  intermediate-scale  cy- 
clones while  Rasmussen  observed  strong  latent  heat  release 
and  a  tropical  structure.   Both  observed  characteristics  are 
present  in  the  diabatic  experiment  and  appear  to  contribute 
to  the  growth  of  the  intermediate-scale  waves. 

B.   SURFACE  FLUX  MODIFICATION  OF  GROWTH 

The  surface  fluxes  during  the  initial  growth  phase  are 
very  small  in  terms  of  normal  atmospheric  values  (Table  4-2). 
This  is  due  to  the  specification  of  the  initial  fields  and 
also  because  the  surface  wind  speed  associated  with  the  small 
amplitude  perturbation  is  still  small.   However,  it  will  be 
shown  that  these  small  surface  fluxes,  through  interaction 
of  the  PBL  with  the  model  diabatic  processes,  have  a  signi- 
ficant effect  on  the  model  atmosphere. 

Table  4-2 

Surface  fluxes  of  zonal  mean  sensible  heat 
(S)  and  moisture  (Q)  in  W/m^  at  day  2. 
(Q  is  in  latent  heat  equivalent.) 

Latitude    Flux     Experiment 

D        D-FS        D-FQ 
46°N        S        18.3      -  25.6 

Q        48.3     42.0 

36°S        S        15.2      -  39.5 

Q        42.5     34.4 

The  static  stability  of  the  surface  flux  experiments 
decreases  in  the  same  sense  as  in  D  (Fig.  4-4  a-b) .   D-FS 
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Figure  4-4.   As  in  Fig.  4-3.   S  and  Q  represent  experi' 
ments  D-FS  and  D-FQ  respectively. 
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(no  surface  sensible  heat  flux)  follows  D  closer  than  does 
D-FQ  (no  surface  moisture  flux) .   Comparison  with  experiments 
A  and  A36  (Fig.  4-3)  indicates  that  static  stability  in  ex- 
periments D-FS  and  D-FQ  follows  that  of  A36  more  closely  than 
that  of  A.   Consequently,  the  reduction  in  static  stability 
should  increase  the  growth  of  the  shorter  waves  in  D-FS  and 
D-FQ  to  a  level  at  least  comparable  to  A36.   This  shift,  how- 
ever, should  occur  gradually  as  static  stability  is  reduced 
and  not  during  the  initial  24  to  48  h  of  integration. 

Figs.  4-5  a-b  give  the  amplitudes  of  the  meridional  wind 
components  at  day  2  and  day  3  by  wavenumber  for  the  diabatic 
experiments.   In  the  preceding  section  wave  growth  rates  be- 
tween days  2.5  and  3.5  are  derived  for  the  adiabatic/diabatic 
comparison.   Here  the  actual  amplitudes  indicate  the  temporal 
change  of  wave  growth  in  response  to  static  stability  and 
other  factors.   Note  also  that  values  given  in  Fig.  4-5  are 
vertically  averaged,  whereas  the  values  in  Fig.  4-1  are  at  a 
single  level.   The  values  for  wave  number  6  are  much  larger 
than  the  other  wave  numbers  at  day  2  because  the  model  was 
initialized  with  a  value  of  approximately  .5  m/s  in  wave  num- 
ber 6.   This  results  in  a  larger  amplitude  for  this  wavenumber 
even  though  it  may  have  a  smaller  growth  rate  than  the  higher 
wave  numbers.   The  exponential  growth  rates,  which  depend  on 
the  ratios  of  the  day  3  to  the  day  2  amplitudes,  are  not 
shown  here. 

The  largest  scale  wave,  wave  number  6,  is  relatively 
insensitive  to  the  removal  of  the  surface  fluxes.   D-FS  (no 
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Figure   4-5. 


Amplitude  (m/s)  of  the  meridional  wind  per- 
turbation by  wavenumber  for  day  2  (stippled) 
and  day  3  (hatched)  for   a.  46°N  and 
b.  36°S.   D,  S  and  Q  as  in  Fig.  4-4. 
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surface  sensible  heat  flux)  and  D-FQ  (no  surface  moisture 
flux)  exhibit  10%  to  15%  less  amplitude  in  wave  number  6  than 
does  D.   There  is  also  a  smaller  growth  rate  in  the  fall  case 
compared  to  the  spring  case  with  the  stronger  jet.   The  re- 
duction in  amplitude  in  D-FS  differs  from  the  study  of  dia- 
batic  effects  on  baroclinic  instability  discussed  in  Chapter 
III.   This  may  be  due  to  the  cooperative  reduction  in  surface 
moisture  flux  in  this  experiment  as  noted  in  the  results  pre- 
sented in  Table  3-2. 

Wave  number  12  is  the  slowest  growing  wave  during  the 
first  one  to  two  days  in  all  the  diabatic  experiments.   Be- 
yond this  time,  however,  it  becomes  the  fastest  growing  wave 
in  experiments  D  and  D-FS.   Referring  to  Fig.  4-4,  the  corre- 
lation with  the  decrease  in  low-level  static  stability  appears 
to  be  very  high. 

Wave  number  18  shows  the  greatest  initial  growth  in  the 
diabatic  experiments.   This  wavelength  and  wavenumber  24  are 
also  the  most  sensitive  to  the  surface  fluxes  of  sensible  heat 
and  moisture.   The  amplitude  of  wavenumber  18  at  days  2  and 
3  in  D-FS  is  20%  to  40%  less  than  in  D.   For  D-FQ,  the  ampli- 
tudes in  both  hemispheres  are  less  than  50%  of  those  in  D. 
Wave  number  18  growth  is  more  rapid  in  the  fall  or  weaker  jet 
hemisphere  in  the  diabatic  experiments.   This  behavior  is 
opposite  to  that  of  wave  number  6,  which  grows  more  rapidly 
in  the  spring  hemisphere.   This  difference  of  growth  could  be 
either  a  response  to  lower  values  of  vertical  wind  shear  as 
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found  by  Rao  and  Ferreira  (1979),  or  a  response  to  higher 
sea-surface  temperatures.   Considering  changes  in  growth 
rate  of  wavenumber  18,  the  greatest  decrease  from  loss  of 
the  surface  sensible  heat  flux  during  the  first  two  days  oc- 
curs in  the  fall  sector.   However,  wavenumber  18  grows  at 
the  same  rate  between  days  2  and  3  in  the  fall  case  for  all 
three  experiments,  which  implies  that  some  other  process  has 
become  dominant. 

The  vertical  distribution  of  the  meridional  velocity  per- 
turbation for  experiments  D,  D-FS,  D-FQ  and  A  is  given  in 
Figs.  4-6  a-b.   A  Fourier  analysis  routine  is  used  to  filter 
out  the  mean  meridional  velocity  and  all  perturbation  com- 
ponents with  a  wavenumber  higher  than  30.   Several  tendencies 
already  noted  in  Figs.  4-1  and  4-5  are  present  here.   The 
largest  amplitude  in  A  is  in  the  spring  hemisphere,  whereas 
it  is  in  the  fall  hemisphere  for  the  diabtaic  experiments. 
The  stronger  spring  jet  controls  the  growth  in  A.   The  pertur- 
bation in  A  is  much  smaller  than  that  found  in  the  diabatic 
experiments — as  much  as  a  factor  of  10  smaller  than  D  at  9  00 
mb  at  4  6°N.   Removal  of  the  surface  fluxes,  D-FS  and  D-FQ, 
makes  the  amplitude  of  the  perturbation  approach  that  of  the 
adiabatic  experiment.   Withholding  the  surface  moisture  flux, 
D-FQ,  has  a  larger  effect  than  removal  of  the  surface  sensible 
heat  flux.   The  contributions  of  the  higher  wave  numbers  (see 
Fig.  4-5),  which  grow  more  rapidly  in  the  fall  hemisphere,  are 
responsible  for  the  larger  amplitude  in  that  hemisphere  of 
the  diabtaic  experiments. 
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Figure  4-6. 


Amplitude  (m/s)  at  day  3  of  the  meridional 
wind  perturbation  summed  over  wavenumbers 
6  through  30  for   a.  46°N  and   b.  36°S. 
Experiments  D,  S(D-FS),  Q(D-FQ)  and  A  are 
compared. 
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The  largest  amplitude  in  D  is  at  the  lowest  model  level. 
This  is  due  to  the  contributions  of  the  high  waven umber  com- 
ponents which  tend  to  be  associated  with  shallow  systems. 
In  D-FQ  and  D-FS ,  the  greater  response  of  the  shorter  waves 
to  the  loss  of  the  surface  fluxes  causes  the  greatest  reduc- 
tion in  growth  at  the  lowest  level.   Again,  D-FS  and  D-FQ  tend 
toward  the  adiabatic  experiment  which  has  its  largest  ampli- 
tude at  30  0  mb. 

The  temperature  perturbations  are  much  larger  in  the 
diabatic  experiments  than  in  A  (Figs.  4-7  a-b) .   The  inclusion 
of  the  surface  fluxes  increases  the  amplitude  of  the  pertur- 
bations.  The  percentage  decrease  in  the  temperature  pertur- 
bation of  D-FS,  D-FQ  and  A  relative  to  D  does  not  reflect  the 
larger  decrease  in  wave  growth  that  is  found  for  the  meridional 
velocity  perturbation.   The  implication  here  is  that  the  con- 
tribution to  growth  of  the  wave  is  not  in  direct  response 
to  the  addition  of  low-level  heat  through  the  surface  fluxes. 
These  fluxes  must  also  indirectly  cause  the  greater  wave 
growth  exhibited  by  the  meridional  velocity  component. 

The  vertical  velocity  perturbations  (Figs.  4-8  a-b,  derived 
as  for  the  meridional  velocity  perturbation)  appear  to  be  the 
link  between  the  temperature  perturbations  and  the  wave  growth 
as  seen  in  the  velocity  perturbations.   The  vertical  velocity 
perturbation  in  D  is  10  to  20  times  that  of  A.   Haltiner 
(1971)  describes  how  larger  vertical  velocities  resulting 
from  diabatic  heating  lead  to  a  greater  conversion  of  eddy 
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available  potential  energy  to  eddy  kinetic  energy,  resulting 
in  enhanced  wave  growth.   A  large  portion  of  the  difference 
between  vertical  velocities  in  D  and  A  is  clearly  due  to  the 
surface  fluxes  of  sensible  heat  and  moisture,  with  surface 
moisture  flux  (D-FQ)  having  the  greatest  effect.   The  hemis- 
pheric difference  is  also  evident  here.   At  46°N  with  a 
higher  SST  and  weaker  jet,  D-FS  has  a  vertical  velocity  per- 
turbation 40%  to  50%  less  than  that  in  D.   In  the  spring 
hemisphere,  the  difference  is  less  than  about  10%.   The  slightly 
larger  value  at  400  mb  in  the  spring  case  becomes  slightly 
less  than  in  D  when  the  values  are  averaged  over  34°  to  36°S. 
In  the  spring  hemisphere,  vertical  velocities  associated  with 
the  stronger  jet  may  be  sufficient  to  trigger  the  release  of 
latent  heat.   In  the  fall  hemisphere  with  the  weaker  jet,  the 
contribution  of  the  larger  surface  sensible  heat  fluxes  in 
generating  vertical  velocities  may  be  a  major  factor  in  trigger- 
ing the  release  of  latent  heat.   Hemispheric  differences  are 
also  evident  in  D-FQ  (no  surface  moisture  flux) .   In  the  fall 
case,  vertical  velocity  perturbations  range  from  an  approxi- 
mately 85%  reduction  from  those  in  D  at  200  mb  to  a  65%  reduc- 
tion at  800  mb.   The  corresponding  spring  case  range  is  70% 
to  one  of  40%.   The  surface  sensible  heat  and  moisture  fluxes 
appear  to  be  cooperative  (as  in  Table  4-2) ,  and  the  lower  SST 
of  the  spring  hemisphere  results  in  a  smaller  moisture  flux 
or  influence  by  the  moisture  flux. 

The  contributions  of  the  surface  fluxes  to  the  distribu- 
tion of  the  vertical  velocity  is  not  the  same  for  both  moisture 
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flux  and  sensible  heat  flux.   In  D-FQ,  the  largest  perturba- 
tion is  in  the  lowest  layer,  800  mb,  in  both  hemispheres. 
This  is  due  to  the  surface  sensible  heat  flux  causing  the 
greatest  heating  at  this  level.   In  D-FS,  the  maximum  is  at 
600  mb,  due  to  heating  caused  by  the  release  of  latent  heat 
at  higher  levels.   Again,  a  cooperation  between  the  two  sur- 
face fluxes  may  be  required  to  produce  the  large  vertical 
velocities  associated  with  D. 

In  summary,  the  surface  fluxes  of  sensible  heat  and  mois- 
ture contribute  to  a  reduction  in  low-level  static  stability 
in  the  diabatic  model  atmosphere.   This  reduction  is  suffi- 
cient to  shift  the  wavelength  of  maximum  growth  in  the  model 
atmosphere  to  shorter  wavelengths  than  those  found  in  the 
adiabatic  model  atmosphere.   In  addition,  latent  heating  in 
the  diabatic  experiments  enhances  the  development  of  shorter 
wavelengths.   This  latent  heating  is  very  sensitive  to  the 
availability  of  low-level  moisture  and,  to  a  lesser  degree, 
surface  sensible  heat  fluxes.   The  fall  case  with  a  weaker 
jet  and  higher  sea-surface  temperatures  is  much  more  sensi- 
tive to  changes  in  the  surface  fluxes. 
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V.   SECONDARY  LOW  DEVELOPMENT  EXPERIMENT 

Several  observational  studies  of  secondary  lows  or  polar 
lows  have  been  performed.   Nitta  and  Yamamoto  (1974)  and 
others  investigated  the  development  of  weak,  small-scale 
cyclones  over  the  East  China  Sea.   Mullen  (1979)  and  Reed 
(1979)  have  completed  case  studies  of  polar  lows  which  form 
in  polar  outbreaks  associated  with  major  frontal  systems  or 
large-scale  cyclones  over  the  eastern  Pacific  Ocean.   Rasmussen 
(1979)  has  investigated  small-scale,  intense  cyclones  in  the 
North  Atlantic  Ocean  which  he  also  terms  polar  lows.   These 
systems  have  two  common  features.   They  are  small  or  inter- 
mediate scale  (zonal  wavelength  of  1000  to  2000  km),  and  they 
originate  over  the  ocean.   Reed  (19  79)  has  also  noted  that  very 
few  intermediate-scale  cyclones  occur  over  land.   Because 
these  lows  appear  mainly  over  the  ocean,  it  is  felt  that  air- 
sea  fluxes  are  likely  to  be  important  during  their  evolution. 

Intermediate-scale  waves  are  present  throughout  the  dia- 
batic  model  experiments,  although  they  are  dominated  by  the 
large-scale  primary  cyclone.   This  chapter  examines  the  evo- 
lution of  the  secondary  cyclones  using  the  technique  described 
above  of  selectively  removing  the  surface  flux  processes  from 
the  diabatic  model.   Table  5-1  is  a  list  of  the  integrations 
performed  for  this  experiment.   Fields  from  the  full  diabatic 
experiment  are  used  to  initialize  the  other  model  experiments 
at  day  4  (S4)  and  day  9  (S9)  . 
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Table  5-1 
Experiments  for  the  secondary  development  study 


Exp. 

Initial 
Fields 

Integration 
Period 

D 

S  (wave  no.  6) 

15  days 

D-FS(4) 

S4  (fields  from  D 
at  day  4) 

6  days 

D-FQ(4) 

S4 

6  days 

D-FSQ(4) 

S4 

6  days 

F(4) 


S4 


6  days 


Model 
Version 

complete  diabatic 

diabatic  w/o  surface 
sensible  heat  flux 

diabatic  w/o  surface 
moisture  flux 

diabatic  w/o  surface 
moisture  or  sensible 
heat  fluxes 

adiabatic  w/friction 


D-FS(9)    S9  (fields  from  D 
at  day  9) 

D-FQ(9)    S9 


D-FSQ(9)   S9 


F(9 


S9 


3  days     diabatic  w/o  surface 
sensible  heat  flux 

3  days     diabatic  w/o  surface 
moisture  flux 

3  days     diabatic  w/o  surface 
moisture  or  sensible 
heat  fluxes 

3  days     adiabatic  w/friction 


See  Sections  II. C  and  II. D  for  a  more  complete  dis- 
cussion of  model  variations  and  initial  conditions. 
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Four  cases  involving  development  of  a  secondary  cyclone 
wave  are  examined.   Surface  pressure  maps  for  an  illustra- 
tive time  from  each  case  are  presented  in  Figs.  5-1  a-d. 
The  first  three  maps  are  from  D  and  the  last  is  from  D-FSQ. 
A  common  feature  of  the  first  three  cases  is  that  the  second- 
ary lows  develop  in  the  equatorward  flow  over  warmer  water, 
similar  to  the  polar  lows  studied  by  Reed  (1979)  and  Rasmussen 
(1979).   The  first  two  cases  examined  (see  Figs.  5-la  and  b) 
are  fall  hemisphere  developments,  while  the  last  two  cases 
(Figs.  5-lc  and  d)  are  spring  hemisphere  developments. 
Intermediate-scale  waves  grow  more  readily  in  the  initial  de- 
velopment experiments.   Not  only  is  the  number  of  secondary 
lows  larger  in  the  fall  hemisphere  than  in  the  spring  hemis- 
phere of  the  diabatic  experiment,  they  are  also  more  intense. 
It  will  be  shown  that  the  secondary  lows  in  the  fall  case  may 
evolve  into  the  next  primary  cyclone.   In  the  spring  case, 
the  secondary  lows  remain  relatively  weak  compared  to  the 
primary  cyclone  and  dissipate  more  rapidly. 

Expanded  analyes  of  various  model  atmosphere  parameters 
and  derived  quantities  are  presented  below.   Analyses  cover 
an  area  of  18°  lat  by  22.5°  long  or  10  x  10  model  grid  points. 
The  expansion  is  about  the  same  center  point  for  any  given 
time  and  hemisphere,  and  is  chosen  as  approximately  the  posi- 
tion of  the  center  of  the  secondary  cyclone  in  D-FS ,  except 
as  noted.   Southern  hemisphere  analyses  are  inverted  with 
south  at  the  top  of  the  map  for  ease  of  comparison  with  the 
northern  hemisphere  cases. 
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30 
LONGITUDE 


Figure  5-1. 


Surface  pressure  maps  for  D  for   a.  Case  A, 
day  6.5,   b.  Case  B,  day  11,   c.  Case  C, 
day  11  and   d.  Case  D,  D-FSQ,  day  10.5. 
Contour  interval  4  mb.   Secondary  low  for 
case  study  indicated  by  an  arrow  in  each 
figure . 
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A.   CASE  A — FALL  HEMISPHERE  DAYS  6-8 

Case  A  gives  an  example  of  selective  development  of  a 
secondary  low.   At  day  6  of  D  a  double  low  exists  approximately 
35  degrees  west  of  the  main  primary  low  (see  position  in  Fig. 
3-2a) .   A  very  similar  system  exists  in  the  same  location  in 
experiments  D-FS,  D-FQ  and  D-FSQ  (see  Figs.  3-7  a-c,  respec- 
tively) .   At  day  6  all  four  diabatic  experiments,  D,  D-FS, 
D-FQ  and  D-FSQ,  exhibit  a  dual  low  structure  with  a  1004  mb 
central  pressure  (see  expanded  surface  pressure  analyses 
Figs.  5-2  a-d) .   These  secondary  lows  evolve  differently  dur- 
ing the  next  12  h  (Figs.  5-3  a-d  and  5-4  a-d).   The  cyclones 
in  D-FS  and  D-FSQ  (no  surface  sensible  heat  flux)  are  more 
intense  than  their  counterparts,  D  and  D-FQ.   In  addition, 
in  D-FQ  and  D-FSQ  (no  surface  moisture  flux)  the  westward 
lobe  of  the  initial  dual  low  (see  Fig.  5-2a)  has  developed, 
causing  a  phase  lag  of  9°  long  behind  the  lows  in  D  and  D-FS 
after  12  h.   A  slight  phase  lag  (1°  long)  is  also  associated 
with  the  addition  of  the  surface  sensible  heat  flux.   The  low 
in  D  lags  that  in  D-FS  and,  similarly,  the  low  in  D-FQ  lags 
that  in  D-FSQ. 

In  D  and  D-FS,  the  secondary  cyclone  intensifies  as  it 
moves  into  the  northerly  flow  region  of  the  primary  cyclone 
(see  Figs.  5-5  a-d),  and  eventually  absorbs  and  replaces  the 
primary  cyclone.   This  is  not  true  in  D-FQ  and  D-FSQ.   In 
the  latter  experiments  the  secondary  cyclone  gradually  weakens 
and  is  absorbed  by  the  primary  cyclone.   The  secondary  lows 
are  relatively  sensitive  to  the  air-sea  fluxes.   Therefore, 
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modification  of  the  fluxes  can  result  in  a  greatly  different 
evolution  in  the  atmosphere  on  a  time  scale  of  three  to  four 
days.   In  this  case,  the  large  differences  are  due  more  to 
selection  of  a  different  perturbation,  rather  than  being  due 
to  significant  differences  in  intensification  rates. 

The  major  differences  in  development  occur  in  response  to 
the  low-level  temperature  advection  and  the  release  of  latent 
heat.   The  low  develops  within  a  strong  baroclinic  zone  (see 
Figs.  5-6  a-d)  with  the  strongest  temperature  gradient  and 
thermal  advection  located  to  the  east  of  the  low.   The  low- 
level  temperature  gradient  is  largest  in  D-FSQ  (no  surface 
flux)  and  smallest  in  D  (full  diabatic) . 

Precipitation  is  restricted  to  the  region  east  of  the 
secondary  lows  in  all  of  the  experiments  (see  Figs.  5-7  a-b 
and  5-8  a-d)  .   At  day  6  (Fig.  5-7)  ,  the  precipitation  rate 
associated  with  the  developing  secondary  low  in  D  and  D-FS  is 
up  to  11  cm/day,  although  this  rate  decreases  slightly  at  day 
6.25  (Fig.  5-8).   The  precipiation  rate  in  D-FQ  and  D-FSQ 
is  less  than  2  cm/day  at  day  6  (not  shown)  and  increases 
slightly  to  2  and  5  cm/day,  respectively,  at  day  6.25.   The 
maximum  precipitation  occurs  over  the  region  of  maximum 
horizontal  temperature  advection  (Fig.  5-6)  in  all  experiments 
Convective  precipitation  is  less  than  2  0%  of  the  total  precipi- 
tation in  all  cases.   The  low-level  specific  humidity  fields 
at  day  6  (Fig.  5-9  a-d)  reflect  the  loss  of  the  surface 
moisture  flux  in  D-FSQ  and  D-FQ.   Specific  humidity  is  approxi- 
mately 30%  higher  in  the  region  of  maximum  precipitation  in 
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Figure    5-7. 


Precipitation  rate  for  day  6  and  expanded 
region  as  in  Fig.  5-2  for   a.  D  and   b.  D-FS 
Contour  interval  2  cm/day.   Experiments  D-FQ 
and  D-FSQ  not  shown  as  rate  less  than  2 
cm/day . 
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D  and  D-FS.   The  pattern  of  maximum  specific  humidity,  and 
also  the  low-level  wind  field  seems  to  indicate  strong  mois- 
ture advection  from  the  southwest  of  the  secondary  low.   In 
D-FQ  and  D-FSQ  only,  a  region  of  relatively  low  specific 
humidity  is  present  in  the  western  quadrant  of  the  secondary 
low. 

The  direction  of  motion  of  the  secondary  low  is  toward 
the  southeast  in  all  of  the  experiments.   This  is  toward  the 
maximum  low-level  thermal  gradient,  the  area  of  highest  precipi- 
tation and  largest  low-level  moisture.   The  implication  is 
that  the  development  is  brought  about  by  low-level  thermal 
processes  as  noted  by  Petterssen  et  a_l .  ,  (1962),  and  that  the 
surface  fluxes  of  sensible  heat  and  moisture  significantly 
modify  this  development. 

To  investigate  the  relative  importance  of  the  surface 
fluxes  in  the  development  of  the  secondary  cyclone,  the  ther- 
mal terms  of  Petterssen 's  development  equation  (Petterssen, 
1956)  are  evaluated.   Petterssen 's  development  equation  de- 
rived in  a  sigma  coordinate  system  is: 


HJ3l    m     Hsu  _  g  v2  3h 


3t      at     f    at  ' 


B 


where: 


5    =   geostrophic  vorticity 
g    =   acceleration  of  gravity 
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f   =   coriolis  parameter,  and 

h  =   thickness  from  level  0  to  level  u. 

Here,  term  A  represents  the  local  change  in  geostrophic 
vorticity  at  the  surface,  term  B  is  the  local  change  in  geo- 
strophic vorticity  at  an  upper  level  and  term  C  is  the  change 
in  thermal  vorticity  due  to  a  change  in  thickness  caused  by 
warming  or  cooling  the  column  of  air. 

Terms  A  and  B  can  be  inferred  from  the  surface  pressure 
fields  and  upper-level  wind  fields.   Term  C  can  be  expanded 
to  evaluate  the  contributions  of  the  individual  thermodynamic 
processes  as  follows: 


E 


+   r~  dT}  6  (ln  a)  '  (5"2) 

p 


where 


Vu   =   horizontal  wind  on  a  sigma  surface 

rl 

V    =   horizontal  gradient  operator  on  a  sigma 
surface 

T   =   temperature 

tt   =   surface  pressure 
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a    =   specific  volume 

C   =   specific  heat  at  constant  pressure 
P 

a    =   p/tt  where  p  is  pressure  on  a  sigma  surface 

a    =   da/dt 

dH/dt  =   heating  rate  due  to  diabatic  effects  such 
as  radiation  and  latent  heat  release,  and 

R   =   gas  constant  for  dry  air. 

Here,  term  D  represents  the  change  in  thickness  due  to  tempera- 
ture advection,  term  E  the  change  due  to  dynamic  warming  or 
cooling  and  term  F  the  change  due  to  diabatic  heating.   In 
general,  warm  advection  and  precipitation  (terms  D  and  F) 
tend  to  enhance  low-level  vorticity,  whereas  the  compensating 
upward  vertical  motion  (term  E)  tends  to  oppose  development. 
Certain  assumptions  are  necessary  to  evaluate  the  terms 
of  the  above  equations.   The  upper  level  is  chosen  as  sigma 
level  4  which  is  approximately  500  mb .   Only  a  very  weak  dis- 
turbance in  essentially  straight  flow  exists  at  this  level. 
The  equations  were  also  evaluated  at  one  of  the  time  periods 
using  sigma  level  5  as  the  upper  level  (approximately  700  mb) 
and  yielded  qualitatively  the  same  results.   The  lowest  model 
layer  contains  an  idealized  planetary  boundary  layer.   Cal- 
culations in  this  layer  are  performed  using  an  average  of 
the  sigma  level  6  (900  mb)  values  and  the  surface  values  (19.5 
m  for  wind  values) ,  rather  than  recalculate  the  complete 
boundary  layer  formulation. 
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Term  B  of  the  Petterssen  development  equation  represents 
the  contribution  to  low-level  vorticity  due  to  the  change  in 
upper-level  vorticity.   Comparison  of  sigma  level  4  (500  mb) 
winds  from  D  and  D-FQ  at  days  5.5,  6  and  6.2  5  revealed  nearly 
identical  (±  1  m/s)  winds  with  a  small  dual  trough  located 
vertically  above  the  surface  secondary  lows  in  both  experi- 
ments.  The  contribution  to  the  secondary  development  from 
term  B  appears  to  be  minor  and  in  any  case  should  be  similar 
in  the  two  experiments. 

Term  D  represents  the  change  in  thickness  due  to  horizon- 
tal temperature  advection.   The  region  of  strongest  advection 
is  discussed  above  in  terms  of  the  low-level  temperature  and 
wind  fields  (Fig.  5-6).   Term  D  (see  Figs.  5-10  a-d)  gives 
the  thickness  change  in  the  lower  model  atmosphere  produced 
by  the  horizontal  temperature  advection.   The  smallest  thick- 
ness increases  east  of  the  secondary  lows  occur  in  D,  and 
the  largest  in  D-FSQ.   The  area  of  positive  thickness  increases 
is  also  much  larger  in  D-FS,  D-FQ  and  D-FSQ  (no  surface  flux). 
This  maximum  thickness  increase  occurs  five  degrees  westward 
in  D-FQ  and  D-FSQ  (no  surface  moisture  flux  experiments) . 
This  is  consistent  with  the  development  of  the  western  per- 
turbation in  these  experiments. 

Term  E  (Figs.  5-11  a-d)  represents  the  dynamic  cooling 
of  the  rising  air  and  acts  as  a  braking  term  to  oppose  develop- 
ment.  This  term  is  influenced  both  by  the  lapse  rate  of  the 
environmental  air  and  by  the  magnitude  of  the  vertical  motion-- 
either  larger  upward  motion  or  a  smaller  lapse  rate  leads  to 
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greater  braking.   D  and  D-FS  have  very  large  values  of  term 
E.   The  secondary  low  in  D  has  the  highest  central  pressure 
at  day  6.25  (see  Fig.  5-3) ,  and  yet  the  magnitude  of  term  E 
is  twice  that  of  D-FQ  and  D-FSQ.   As  discussed  in  the  initial 
growth  experiments,  this  is  associated  with  vertical  veloci- 
ties that  are  two  to  three  times  larger  in  D  and  D-FS. 

Term  F  (Fig.  5-12  a-d)  represents  the  warming/cooling  of 
the  column  of  air  through  the  diabatic  processes  of  radiation, 
precipitation  (latent  heat  release)  and  surface  sensible  heat 
flux.   The  greater  warming  in  D-FS  relative  to  D  seems  to  be 
related  to  the  20%  greater  precipitation  (see  Fig.  5-7)  in 
D-FS  at  day  6.25.   This  difference  obscures  the  influence  of 
the  surface  sensible  heat  flux  between  D  and  D-FS. 

Diabatic  heating  (term  F)  contributes  more  to  the  develop- 
ment of  the  secondary  low  than  does  thermal  advection  (term 
D) .   However  the  mechanism  for  lifting  the  air  and  causing 
the  diabatic  heating  is  clearly  the  thermal  advection.   This 
is  because  of  the  near  colocation  of  the  maxima  in  these  fields 
in  Fig.  5-10  and  Fig.  5-12  (also  see  Figs.  5-6  and  5-8).   The 
loss  of  the  surface  fluxes  of  sensible  heat  and  moisture  lead 
to  a  larger  contribution  from  thermal  advection.   The  loss  of 
the  surface  moisture  flux,  however,  has  a  much  greater  effect, 
because  it  influences  both  the  location  of  phasing  as  well  as 
the  intensity  of  the  development.   The  total  contribution  of 
the  thermal  processes  (term  C)  to  the  thermal  vorticity  (Fig. 
5-13)  is  several  times  larger  in  the  experiments  with  surface 
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moisture  flux.   This  process  appears  to  be  the  reason  why 
the  secondary  lows  continue  to  develop  in  D  and  D-FS  while 
they  slowly  dissipate  in  D-FQ  and  D-FSQ. 

Although  there  is  some  deep  convection,  precipitation 
caused  by  large-scale  lifting  appears  to  be  responsible  for 
the  development  of  the  secondary  lows  in  D  and  D-FS.   This 
lifting  is  caused  by  warm  advection  and  occurs  in  a  region 
of  strong  moisture  gradient.   This  moisture  gradient  is  main- 
tained by  the  advection  of  moisture  and  by  the  surface  mois- 
ture flux  and  results  in  much  greater  precipitation  in  D  and 
D-FS. 

B.   CASE  B — FALL  HEMISPHERE  DAYS  10-12 

The  secondary  development  in  this  case  is  similar  to  the 
one  examined  in  Case  A.   Model  versions  D-FS,  D-FQ  and  D-FSQ 
are  also  initialized  with  fields  from  D  at  day  9.   A  1008  mb 
secondary  low  first  appears  at  approximately  day  10  in  all 
model  versions.   At  day  10.5,  1.5  days  after  initialization, 
the  secondary  low  enters  the  southwest  quadrant  of  the  primary 
low  (see  Fig.  5-14  a-d  and  also  Fig.  5-lb)  with  approximately 
the  same  central  pressure,  1007  mb  ±  1  mb ,  and  location  in 
all  the  experiments.   During  the  next  18  h,  the  secondary  low 
develops  (see  Figs.  5-15  a-d  and  5-16  a-d)  much  more  rapidly 
in  D  and  D-FS,  and  becomes  the  primary  low  by  day  12.   In 
D-FQ  and  D-FSQ  (no  surface  moisture  flux)  the  secondary  low 
slowly  dissipates  after  day  11.   Also  the  phase  speed  of  the 
secondary  low  is  less  in  D-FQ  and  D-FSQ.   This  results  in  a 
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phase  lag  of  6°  long  relative  to  the  secondary  low  in  D  during 
the  period  from  day  10.5  to  day  11.25. 

Case  B  closely  parallels  Case  A.   To  simplify  the  dis- 
cussion, only  the  contributions  of  terms  D,  E,  and  F  in  the 
Petterssen  development  equation  will  be  presented.   Term  D 
(Fig.  5-17  a-d)  yields  similar  contributions  as  in  Case  A. 
D-FS  and  D-FSQ  (no  surface  sensible  heat  flux) ,  have  larger 
contributions  to  development  than  the  experiments  with  sensible 
heat  flux,  D  and  D-FQ,  respectively.   Term  E  (Fig.  5-18  a-d) 
again  represents  a  much  larger  braking  in  D  and  D-FS  due  to 
the  larger  vertical  motions  associated  with  the  secondary  low. 
The  diabatic  heating  term  (term  F,  Fig.  5-19  a-d),  represents 
the  largest  contribution  to  the  thickness  increase.   The  sum 
of  terms  D  and  F  is  not  offset  by  term  E  in  D  and  D-FS,  which 
allows  the  secondary  low  to  develop  in  D  and  D-FS,  but  dissi- 
pate in  D-FQ  and  D-FSQ.   The  thermal  advection  term  contributes 
to  larger  development  in  D-FS  relative  to  D  and  less  dissi- 
pation in  D-FSQ  than  in  D-FQ.   Lifting  due  to  warm  advection 
is  clearly  the  dominant  mechanism  in  producing  precipitation, 
and  thermal  advection  is  promoted  by  the  absence  of  the  sur- 
face sensible  heat  flux. 

C.   CASE  C — SPRING  HEMISPHERE  DAYS  10-12 

Case  C  is  an  example  of  a  secondary  low  that  fails  to 
deepen.   As  in  Case  B,  model  versions  D-FS,  D-FQ  and  D-FSQ 
are  initialized  with  fields  from  D  at  day  9 .   A  secondary  low 
is  evident  in  all  model  versions  about  day  10.5,  or  1.5  days 
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after  initialization.   The  low  is  in  approximately  the  same 
location  (see  Fig.  5-20  a-d  and  Fig.  5-lc)  with  only  small 
differences  in  initial  intensity.   The  secondary  low  moves 
through  the  region  southwest  (would  be  northeast  if  plotted 
normally)  of  the  primary  low,  which  corresponds  to  the  loca- 
tions of  the  developments  that  occurred  in  Cases  A  and  B.   In 
this  case,  however,  the  low  merely  maintains  the  same  inten- 
sity for  approximately  12  h,  and  then  slowly  dissipates  in 
all  model  versions  (Figs.  5-21  and  5-22) . 

The  Petterssen  development  equation  terms  have  been  evalu- 
ated for  Case  C  at  day  10.75  (Figs.  5-23,  5-24  and  5-25). 
Term  D  (Fig.  5-23)  reflects  a  somewhat  weaker  thermal  advec- 
tion  contribution  than  in  Cases  A  and  B,  although  the  decrease 
in  Case  C  is  less  than  10%  for  D  and  D-FS.   It  is  the  second- 
ary lows  in  these  two  model  versions  which  develop  strongly 
in  Cases  A  and  B.   Term  F  (Fig.  5-25)  demonstrates  a  decrease 
of  60%  to  80%  in  the  diabatic  warming  of  Case  C  from  Cases  A 
and  B  (Figs.  5-12  and  5-19)  in  experiments  D  and  D-FS.   The 
maximum  values  of  the  diabatic  heating  in  Case  C  are  only  of 
the  same  order  as  those  in  D-FQ  and  D-FSQ  of  Cases  A  and  B. 
The  secondary  low  in  D-FQ  and  D-FSQ  dissipated  eventually  in 
both  Cases  A  and  3.   The  implication  here  is  that  a  large 
contribution  from  diabatic  heating,  term  F,  is  required  for 
the  secondary  low  to  continue  to  develop. 

D.   CASE  D — SPRING  HEMISPHERE  DAYS  9-11 

This  case  differs  from  the  previous  three  cases  in  that 
the  secondary  low  develops  in  the  poleward  flow  in  the  eastern 
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quadrant  of  the  primary  cyclone  (see  Fig.  5-ld) .   Also,  in 
apparent  contradiction  to  the  first  three  cases  and  the  re- 
sults discussed  in  Chapter  III,  the  largest  development  occurs 
in  the  adiabatic  experiment,  F,  with  the  least  development 
in  D. 

Figs.  5-26  a-b  present  the  expanded  surface  pressure  maps 
for  F  corresponding  to  those  presented  in  Figs.  5-20  a-d  and 
5-22  a-d  for  the  diabatic  experiments.   The  secondary  low 
reaches  its  maximum  intensity  between  days  10.5  and  11  in 
all  experiments,  prior  to  entering  the  region  of  strong 
equatorward  flow  where  development  occurs  in  Cases  A  and  B. 
At  this  time,  day  10.5,  the  secondary  low  is  deepest  in  F, 
followed  by  D-FSQ,  and  with  the  highest  central  pressure 
occurring  in  D. 

Temperature  fields  for  experiments  D  and  F  (see  Figs. 
5-27  a-b)  reveal  a  much  larger  low-level  temperature  gradi- 
ent to  the  east  of  the  primary  cyclone  in  F  than  D.   The 
gradients  for  D-FSQ,  D-FS  and  D-FQ  (not  shown) ,  respectively, 
are  also  larger  than  D  and  the  secondary  lows  are  corres- 
pondingly more  intense  (see  Figs.  5-20  a-d).   Expanded  plots 
of  the  low-level  temperature  fields  (Figs.  5-28  a-b)  reveal 
a  temperature  gradient  in  F  approximately  twice  that  in  D. 
Thus,  the  indication  is  that  low-level  temperature  advection 
is  responsible  for  the  difference  in  development  in  this  case. 

In  terms  of  the  Petterssen  development  equation,  the 
secondary  low  developing  in  F  (adiabatic)  has  no  contribu- 
tion from  diabatic  heating  (term  F) ,  yet  it  develops  more 


118 


LONGITUDE 


LONGITUDE 


Figure    5-26. 


Expanded  surface  pressure  maps  of  secondary 
low  in  F  in  the  spring  hemisphere  for   a.  day 
10.5  and   b.  day  11.   Region  as  in  Figs.  5-20 
and  5-22,  respectively.   Contour  interval 
2  mb. 
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Figure  5-27. 


Southern  hemisphere  temperature  fields  at 
sigma  level  6  (900  mb)  at  day  10.5  for 
a.  D  and  b.  F.   Contour  interval  5°C. 
Location  of  secondary  low  indicated  by  X. 
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Figure    5-28. 


As  in  Fig.  5-27  except  for  expanded  region 
as  in  Fig.  5-26a.   Contour  interval  2°C. 
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Figure    5-29. 


Term  D  of  the  Petterssen  development  equa- 
tion for  day  10.5  and  an  expanded  region 
as  in  Fig.  5-26a  for   a.  D  and   b.  F. 
Contour  interval  5  x  10~5  °c/s. 
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than  the  secondary  low  in  D.   Figs.  5-29  a-b  give  the  con- 
tribution of  the  low-level  temperature  advection  (term  D) 
to  development  of  the  secondary  low  at  day  10.5.   The  region 
of  positive  thickness  change  (negative  values  in  figure)  is 
larger  as  well  as  more  intense  in  F. 

Developments  of  this  type  (Case  D)  in  the  eastern  quadrant 
of  the  primary  cyclone  occur  in  both  hemispheres,  although 
the  intensity  changes  are  generally  not  as  marked  as  in  this 
example.   These  situations,  in  which  the  air-sea  fluxes  and 
diabatic  effects  appear  to  hinder  the  cyclone  development, 
are  more  common  and  more  marked  in  the  coarse-mesh  model 
experiments.   It  appears  that  the  relative  role  of  the  physi- 
cal processes  is  changed  in  the  different  resolutions.   In 
the  northern  or  fall  hemisphere,  the  larger  diabatic  effects 
in  experiments  D  and  D-FS  tend  to  mask  the  contribution  of 
the  stronger  low-level  temperature  advection. 

E .   SUMMARY 

These  cases  clearly  demonstrate  the  sensitivity  of  the 
model  development  of  secondary  lows  to  both  the  surface  fluxes 
of  sensible  heat  and  moisture  and  to  the  baroclinicity  of  the 
model  atmosphere.   The  loss  of  the  surface  sensible  heat  flux 
contributes  to  a  greater  low-level  temperature  gradient 
in  both  the  spring  and  fall  hemispheres.   This  contributes  to 
an  increase  in  low-level  temperature  advection,  and  leads  to 
a  greater  development  in  the  experiments  with  no  sensible  heat 
flux.   The  surface  moisture  flux  is  responsible  for  much 
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larger  differences  in  the  development  of  secondary  lows. 
Here  the  physical  mechanism  is  the  latent  heat  release 
triggered  by  vertical  velocities  induced  by  the  low-level 
temperature  advection. 

In  the  fall  hemisphere  with  the  weaker  jet,  the  secondary 
low  develops  rapidly  in  the  experiments  with  a  surface  mois- 
ture flux,  and  eventually  replaces  the  primary,  large-scale 
cyclone.   In  the  experiments  without  surface  moisture  flux, 
the  secondary  lows  slowly  dissipate.   In  addition,  there  is 
a  six-to-ten  degree  phase  lag  of  the  secondary  lows  in  these 
experiments  relative  to  those  which  had  surface  moisture 
fluxes.   This  phase  lag  occurs  primarily  during  the  12  to 
18  h  period  in  which  the  secondary  lows  develop  most  rapidly. 

In  sharp  contrast,  in  the  spring  hemisphere  with  the 
stronger  jet  and  lower  sea-surface  temperatures,  the  second- 
ary lows  dissipated  after  one  or  two  days,  without  intensi- 
fying greatly.   The  low-level  temperature  advection  in  this 
hemisphere  is  only  about  10%  less  than  in  the  experiments  in 
which  the  secondary  lows  develop  strongly  in  the  fall  hemis- 
phere.  The  contribution  from  the  diabatic  heating,  however, 
is  much  less  in  the  spring  hemisphere  experiments  with  sur- 
face moisture  flux,  and  is  comparable  to  the  experiments  in 
which  the  secondary  low  dissipated  in  Cases  A  and  B. 

Reed  (1979)  concluded  that  the  polar  lows  in  the  eastern 
North  Pacific  Ocean  were  predominately  a  baroclinic  phenomena 
with  possible  contributions  from  thermal  and  barotropic  proc- 
esses.  Rasmussen  (1979),  in  a  similar  study  of  eastern  North 
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Atlantic  Ocean  polar  lows,  concluded  that  they  were  pre- 
dominately a  thermal  instability  phenomena.   The  results  of 
the  model  experiments  presented  here  demonstrate  that  both  of 
the  above  situations  could  exist  in  the  marine  atmosphere. 
The  instability  process  which  dominates  in  any  given  instance 
is  sensitive  to  the  surface  fluxes  of  sensible  heat  and  mois- 
ture, as  well  as  to  the  degree  of  baroclinic  instability 
present  in  the  atmosphere.   In  these  experiments,  the  surface 
fluxes  are  critically  important  on  a  time  scale  of  two  to 
four  days,  which  is  well  within  the  predictability  of  the 
real  atmosphere. 
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VI.   CONCLUSIONS 

The  objective  of  this  research  is  to  use  a  numerical  model 
to  examine  the  role  of  the  air-sea  fluxes  in  extratropical 
cyclogenesis  over  the  open  ocean.   The  method  employed  is  to 
specify  analytic  initial  atmosphere  and  ocean  conditions  which 
are  typical  of  fall  and  spring  conditions  in  the  real  atmos- 
phere.  A  state-of-the-art  atmospheric  prediction  model  is 
used  to  examine  the  contributions  of  the  surface  fluxes  of 
moisture  and  sensible  heat.   This  is  accomplished  by  selec- 
tively removing  the  surface  flux  parameterizations  from  the 
model  and  comparing  the  results  with  a  control  run  generated 
with  the  complete  diabatic  model.   Consequently,  this  re- 
search emphasizes  the  role  of  the  individual  processes  in  the 
model  and  in  the  atmosphere  rather  than  the  response  of  the 
model  to  small  changes  in  the  sea-surface  temperatures  (e.g., 
Spar  and  Atlas,  1975)  . 

The  surface  fluxes  make  a  significant  contribution  in 
these  simulations  of  cyclone  evolution.   In  the  fall  hemi- 
sphere with  a  weaker  upper-level  jet  and  higher  sea-surface 
temperature,  the  removal  of  the  surface  fluxes  results  in  a 
large  reduction  in  the  initial  growth  rate  of  the  model  cy- 
clones.  In  the  presence  of  the  more  intense  jet  and  lower 
sea-surface  temperature  of  the  spring  hemisphere,  removal  of 
the  surface  fluxes  has  a  smaller  impact  on  the  model  cyclone 
growth.   The  contribution  of  the  surface  fluxes  to  the  growth 
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of  a  pre-existing,  large-scale  cyclone  is  found  to  be  rela- 
tively small  as  expected  from  previous  studies. 

The  smaller-scale  secondary  cyclones  respond  rapidly  to 
the  removal  of  the  model  surface  fluxes.   Removal  of  the  sur- 
face moisture  flux  prevents  the  development  of  secondary 
lows  which  might  otherwise  eventually  replace  the  primary 
cyclone,  as  occurs  in  the  complete  diabatic  model.   Removal 
of  the  surface  sensible  heat  flux  results  in  the  development 
of  more  intense  cyclones  although  the  overall  evolution  of 
the  cyclone  is  not  altered  from  that  in  the  complete  dia- 
batic experiment.   The  result  that  different  perturbations 
eventually  become  the  primary  cyclone  after  three  to  four  days 
is  significant.   One  may  conclude  that  extending  the  useful- 
ness of  numerical  model  forecasts  of  open-ocean  cyclones 
beyond  the  present  limit  of  two  to  three  days  will  require 
accurate  specification  of  the  air-sea  fluxes. 

The  experimental  technique  of  removing  the  surface-flux 
parameterizations  from  the  model  has  no  real  world  analog. 
The  modelled  results  can  not  be  verified  with  observations 
and  indicate  only  the  sensitivity  of  the  model  cyclones  to  the 
air-sea  fluxes.   The  ability  to  perform  experiments  within  a 
controlled,  but  admittedly  idealized,  environment  can  provide 
insights,  however,  that  would  otherwise  require  very  large, 
detailed  statistical  studies.   These  results  do  suggest  why 
some  disturbances  over  the  ocean  develop  while  other,  seem- 
ingly similar,  disturbances  do  not  develop. 
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The  surface  fluxes  of  moisture  and  sensible  heat  act  to 
reduce  the  large  gradients  of  low-level  temperature  and  mois- 
ture in  the  model  atmosphere.   These  gradients  occur  in  re- 
gions of  strong  equatorward  (poleward)  advection  associated 
with  the  primary  cyclone,  especially  over  the  zone  of  largest 
sea-surface  temperature  gradient.  Large-scale  lifting  produced 
by  the  low-level  temperature  advection  triggers  the  release 
of  latent  heat  as  hypothesized  by  Petterssen  e_t  al.  ,  (1962)  . 
These  two  processes  jointly  lead  to  a  large  increase  in  ver- 
tical motion  and,  therefore,  cyclone  development.   Removal  of 
the  surface  sensible  heat  flux  allows  a  larger  low-level 
temperature  gradient,  resulting  in  enhanced  vertical  motions 
and  in  increased  development  of  the  cyclones.   Removal  of 
the  surface  moisture  flux  results  in  a  large  decrease  in  pre- 
cipitation or  latent  heat  release  and,  therefore,  a  large 
decrease  in  cyclone  development. 

Previous  studies  have  focused  on  the  region  behind,  or 
westward,  of  the  cyclone,  where  the  strongest  fluxes  occur. 
The  cooperative  effect  of  lifting  produced  by  low-level  tempera- 
ture advection  and  the  release  of  latent  heat,  however,  results 
in  the  largest  cyclone  response  occurring  ahead,  or  eastward, 
of  the  cyclone.   The  results  of  these  experiments  indicate 
that  the  magnitude  of  the  surface  fluxes  in  terms  of  the 
amount  of  energy  transferred  across  the  air-sea  interface  is 
not  as  important  as  the  location  in  which  they  occur.   This 
resolves  the  apparent  contradiction  found  by  Wei  (1979)  and 
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also  implies  that  accurate  specification  of  the  sea-surface 
temperature  distribution  is  crucial  to  the  evolution  of  open- 
ocean  cyclones,  at  least  on  a  medium  range  time  scale  of 
three  to  five  days. 

This  experiment  addresses  open-ocean  conditions  with 
fixed  sea-surface  temperatures  and  moderate  values  of  surface 
moisture  and  sensible  heat  fluxes.   A  variety  of  cyclone  evo- 
lutions similar  to  those  observed  by  Nitta  and  Yamamoto  (1974), 
Reed  (1979)  and  others  is  produced  in  the  model  integrations. 
Much  more  intense  over-ocean  cyclogenesis  has  been  observed 
(Sanders  and  Gyakum,  1980;  Pyke,  1965)  in  association  with 
land-sea  and  ice-sea  boundaries  where  the  atmosphere  is  far 
from  equilibrium  with  the  ocean  surface  fluxes.   Further  re- 
search is  necessary  to  investigate  more  fully  the  role  of 
the  surface  fluxes  in  these  cases  of  more  intense  cyclogene- 
sis.  Numerical  experiments  of  the  type  described  here  can 
be  helpful  in  understanding  the  physical  processes  involved 
in  these  cases,  and  will  hopefully  lead  to  improved  prediction 
of  extratropical  cyclones  over  the  oceans. 
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